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CHAPTER  1 

AN  INTRODUCTION  TO  PREVENTIVE 
MAINTENANCE  POLICIES 

1.1  INTRODUCTION 

During  the  past  two  decades,  there  has  been  a  continuing  interest 
in  studying  maintenance  models  for  systems  with  stochastic  failures  and 
in  their  military,  industrial,  environmental,  and  ecological  applications. 
These  applications  may  include  the  maintenance  of  complex  electronic 
and/or  mechanical  equipment  such  as  computers  and  airplanes,  control 
of  pollutants  in  the  environment,  maintenance  of  the  human  body,  and 
maintenance  of  ecological  balance  in  populations  of  plants  and  animals. 

Barlow  and  Proschan  (1965),  McCall  (1965),  and  Pierskalla  and 
Voelker  (1973)  have  researched  and  surveyed  various  maintenance 
models  in  the  area  of  maintainability. 

A  broad  class  of  what  in  the  area  of  maintainability  is  known  as 
optimal  preventive  maintenance  models  has  received  the  most  attention  in 
the  literature.  An  equipment  whose  failure  rate  increases  with  time 
(e.g.  Weibull,  truncated  normal,  gamma  distributions)  may  break  down 
during  its  actual  operation  more  frequently  as  it  ages.  These  break¬ 
downs  can  often  be  costly  and  a  preventive  maintenance  strategy  can 
cut  down  on  the  irregularity  of  breakdowns  by  carrying  out  regularly 
planned  repairs  before  the  system  fails.  A  preventive  maintenance 
policy  can  also  determine  the  replacement  time  after  which  the  existing 
system  is  no  longer  economical  to  keep.  To  judge  the  quality  of  main- 
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tenance  schedules,  the  decision  can  be  based  on  several  economic 
criteria,  including  minimizing  the  long-term  average  cost  per  unit  time 
and  maximizing  the  steady  state  availability  or  limiting  efficiency  which 
is  the  fraction  of  up-time  over  a  long  time  interval. 

A  substantial  body  of  literature  on  the  preventive  maintenance  of 
stochastically  failing  systems  has  been  devoted  to  the  analysis  of  two 
strategies.  One  is  called  "Age  Replacement"  or  policy  I,  and  the  other 
is  known  as  "Periodic  Replacement  with  Minimal  Repair  at  Failure"  or 
policy  II.  It  is  usually  assumed  that  the  failure  rate  is  a  strictly 
increasing  function  of  time.  In  earlier  works,  generally,  a  preventive 
maintenance  action  refers  to  a  renewal  where  the  system  becomes  as 
good  as  new.  Another  assumption  that  seems  to  be  made  in  the 
literature  most  frequently  is  the  notion  of  zero  down-time  duration  for 
either  repair  or  replacement. 

Here,  we  present  a  brief  discussion  and  literature  review  on 
Policies  I  and  II. 

1.2  POLICY  I 

This  strategy  is  most  useful  in  maintaining  less  complex  equipment, 
such  as  car  engines.  This  strategy  is  defined  as  follows: 

Perform  preventive  maintenance  at  the  time  of  failure  or  if  T 
hours  of  actual  operation  have  elapsed.  Preventive  mainte¬ 
nance  is  then  rescheduled.  This  policy  is  also  referred  to  as 
"Age  Replacement  Policy." 

Barlow  and  Hunter  (1960)  have  studied  this  preventive  maintenance 
policy.  They  discussed  the  existence  and  uniqueness  of  the  optimum 
cycle  T  ,  which  is  the  time  between  two  successive  scheduled  mainte- 
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nance  actions.  Their  objective  is  to  maximize  the  limiting  efficiency 

* 

(steady  state  availability)  if  the  unit  was  replace  at  failure  or  at  T  , 
whichever  occurred  first.  They  have  also  shown  that  the  optimum 
maintenance  time  is  the  solution  to  an  integral  equation. 

Barlow  and  Proschan  (1965)  have  an  equivalent  model  but  their 
criterion  is  to  minimize  the  long-term  expected  average  cost  per  unit 
time,  which  is  the  expected  total  cost  for  a  cycle  divided  by  the  length 
of  that  cycle. 

Glasser  (1967)  has  obtained  solutions  to  the  age  replacement 
problem  for  Weibull,  gamma,  and  truncated  normal  distributions. 

Fox  (1966)  has  used  a  total  discounted  cost  criterion  to  optimize 

* 

the  type  I  policy.  He  has  shown  that  the  optimum  age  of  the  unit  T 
is  a  solution  to  an  integral  equation. 

Most  age  replacement  policies  presented  in  the  literature  consider 
two  costs  for  replacing  failed  and  non-failed  units.  Scheaffer  (1971) 
extends  the  standard  age  replacement  model  by  including  an  age- 
dependent  maintenance  cost.  This  cost  may  reflect  the  increase  in  total 
cost  due  to  depreciation,  wear,  age,  diminishing  productivity  or  reduc¬ 
tion  in  salvage  value  of  the  equipment.  He  proposed  two  cost  functions 
for  which  the  failure  rate  may  also  be  constant  in  order  to  find  an 
optimum  policy  I. 

In  practice,  the  maintenance  is  not  necessarily  the  replacement  of 
the  whole  system,  but  is  often  the  repair  or  replacement  of  part  of  the 
system.  This  kind  of  maintenance  action  does  not  always  renew  the 
system  and  the  mean  lifetime  of  the  repaired  system  is  usually  less  than 
the  new  one.  Tahara  and  Nishida  (1973)  proposed  models  in  which  the 
system  cannot  recover  completely  after  each  repair.  This  notion  has 


been  used  by  Nguyen  and  Murthy  (1981).  They  have  studied  policy  I 
and  policy  II  with  a  basic  assumption  that,  after  each  repair,  the  life¬ 
time  distribution  of  the  system  changes  in  such  a  way  that  its  failure 
rate  becomes  an  increasing  function  of  the  number  of  previous  repairs. 

1.3  POLICY  II 

This  strategy  is  applied  to  very  expensive  and  complex  equipment 
such  as  computers  and  airplanes,  and  it  is  defined  as  follows: 

Perform  preventive  maintenance  on  the  system  after  it  has 
been  operating  a  total  of  T  hours  regardless  of  the  number  of 
intervening  failures.  Preventive  maintenance  is  then  re¬ 
scheduled.  When  a  failure  occurs,  it  is  corrected  by  a 
minimal  repair. 

This  strategy  is  also  known  as  "Periodic  Replacement  with  Minimal 
Repair  at  Failure."  The  action  of  restoring  a  failed  system  to  operation 
without  affecting  its  failure  rate  is  called  minimal  repair.  The  system  is 
usually  completely  renewed  after  a  fixed  number  of  the  periodic  main¬ 
tenance  cycles. 

Most  preventive  maintenance  policies  presented  in  the  literature 
have  studied  the  type  II  policy  together  with  policy  I,  including  Barlow 
and  Hunter  (1961)  and  Barlow  and  Proschan  (1965). 

Bellman  (1865)  and  Descamp  (1965)  applied  dynamic  programming  to 
policy  II  described  above. 

It  is  assumed,  generally,  that  minimal  repairs  performed,  when 
policy  II  is  followed,  have  zero  time  duration.  But  Sivazlian  (1973)  has 
generalized  this  model  by  permitting  a  random  down-time  duration  for 
each  minima]  repair.  This  random  time  has  an  arbitrary  distribution. 
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His  criterion  is  to  minimize  the  long-term  expected  total  cost.  He  also 
derives  the  necessary  and  sufficient  conditions  for  which  this  policy  II 
is  optimal. 

Makabe  and  Morimura  (1963a,  1963b,  and  1965)  and  Morimura  (1970) 
have  studied  another  similar  policy  (policy  III)  along  with  policies  I  and 
II.  Under  policy  III,  the  system  is  replaced  at  the  failure.  The 
intervening  failures  are  corrected  by  minimal  repairs .  They  have 
shown  the  existence  of  the  optimal  policy  III  both  with  respect  to 
limiting  efficiency  and  another  criterion  which  they  called  the  "ft  i- 
tenance  Cost  Rate."  It  is  defined  as  follows: 

[Cost  per  unit  down  time]  x  [expected  fraction  of  down-time] 

+  [expected  cost  of  all  repairs  and  replacements  during  a  unit  time] 

Morimura  (1970)  has  studied  a  more  general  policy  (IV).  Under 

this  policy,  preventive  maintenance  is  performed  when  at  the  time  of 

failure  either  the  total  operating  time  has  exceeded  t  or  when  Ku 

failure  occurs.  Other  failures  are  corrected  with  minimal  repairs.  Note 

that  when  K  approaches  infinity,  the  preventive  maintenance  is  per- 
* 

formed  at  t  and  consequently  we  obtain  the  type  II  policy.  When 
t  =0,  then  the  preventive  maintenance  is  carried  out  at  the  K 
failure  and  we  lead  to  the  type  III  policy. 

As  mentioned  earlier,  the  works  presented  by  Tahara  and  Nishida 
(1973)  and  Nguyen  and  Murthy  (1981)  are  of  great  importance  in 
studying  policies  I  and  II.  The  former  introduces  models  in  which  the 
system  does  not  recover  completely  after  each  repair,  and  the  latter 
relates  the  repaired  system's  failure  rate  to  the  number  of  repairs. 
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1.4  AN  OVERVIEW 

We  consider  optimization  of  schedules  for  maintenance  or  repairs 
when  two  classes  of  maintenance  strategies,  namely  I  and  II,  are  of 
interest.  These  policies  are  said  to  be  optimum  if  they  minimize  the 
expected  long-term  average  cost  per  unit  time. 

Under  policy  I,  repairs  are  performed  at  failure  or  at  Tj  units  of 

qt 

time,  measured  from  the  (i-1)  maintenance  operation,  and  replacement 
occurs  at  t^  defined  by 

N 

tM  =  I  min  (Ti#  x.)  (1.1) 

i=l  1  1 

where  Xj  is  a  random  variable  representing  the  failure  time  from  the 
most  recent  maintenance  action.  There  is  only  one  kind  of  repair 
associated  with  this  policy,  by  which  the  system  improves  but  does  not 
recover  completely.  Our  basic  assumption  is  that  the  lack  of  complete 
recovery  is  due  to  the  system's  previously  expended  lifetime. 

Under  policy  II,  there  are  two  kinds  of  repairs: 

1.  Minor  repairs  to  restore  a  failed  system  to  operation  without 
affecting  its  failure  rate. 

2.  Major  repairs  to  reduce  the  future  number  of  failures. 

Obviously,  in  this  model  the  only  repair  that  may  be  planned  is 
the  major  one  by  which  the  failure  rate  decreases  but  the  system  will 
not  become  as  good  as  new.  The  reason  that  the  regeneration  does  not 
occur  could  be  based  on  the  previous  number  of  repairs,  age,  ineffi¬ 
ciency  of  repairs,  and  many  other  factors.  However,  our  main 
assumption  is  that  the  system's  age  plays  a  significant  role  on  the 
recovery  after  a  major  repair.  Here,  major  repairs  are  scheduled  at 
Ti,T2 . T^  and  replacement  at 
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lN  = 


N 

i=iTN 


(1.2) 


The  novelty  here  is  that  the  failure  rate  after  a  maintenance 
operation  is  a  function  of  the  system's  previously  expended  lifetime. 
This  generalizes  earlier  work  by  Nguyen  and  Murthy  (1981)  on  the 
simpler  case  where  the  future  failure  rate  depends  on  number  of 
previous  repairs,  but  not  on  the  times  when  they  took  place. 

Our  goal  is  to  optimize  both  policies  I  and  II,  separately,  by 
finding  a  set  of  successive  time  intervals  TlfT2, . . .  ,T^  and  the  number 
N,  where  a  replacement  is  made  at  tN  given  by  relations  (1.1)  and  (1.2). 

These  models  are  only  applicable  when  the  system's  failure  rate  is 
increasing  with  time.  All  of  our  results  are  valid  when  the  underlying 
distribution  is  Weibull,  whereas  most  can  be  applied  to  a  broad  class  of 
increasing  hazard  rate  equipments.  The  failure  characteristic  of  the 
system  at  the  beginning  of  the  ith  time  interval  is  described  by  a 
conditional  distribution  function  ^.(tie.^),  where  0.^  describes  the 
previous  history  of  the  system.  Our  contribution  can  be  characterized 
by  letting  0. represent  the  age  factor  affecting  the  system's  failure 
rate.  Policies  I  and  II  are  studied  in  Chapters  2  and  3  respectively. 

Section  2.2  introduces  policy  I,  which  is  a  generalization  of  "repair 
at  failure  or  at  at  time  T,  whichever  comes  first"  to  reduce  the  irregu¬ 
larity  of  breakdowns,  assumptions  involving  zero  time  duration  for  a 
repair  or  replacement,  and  constant  replacement.  Repair  and  break¬ 
down  costs  are  also  discussed  in  this  section. 

Since  we  assume  that  the  recovery  after  each  repair  is  not  fully 
achieved,  function  0-  is  defined  in  Section  2.3  to  relate  the  system's 
age  and  the  number  of  repairs  to  the  failure  rate.  A  formulation  of  the 
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long-term  expected  average  cost  per  unit  time  for  policy  I  is  also  found 
in  this  section. 

Section  2.4  contains  a  general  foundation  of  the  type  I  policy  as  a 
nonlinear  programming  problem  (NPI)  and  a  discussion  on  three  special 
cases : 

1.  Replacement  at  time  T  or  failure,  whichever  occurs  first. 

2.  The  system  ages  rapidly  or  repairs  are  destructive  (91»0O). 

Vei-i 

3.  The  system  ages  slowly  (  -x -  is  small) 

i-I 

A  reformulation  which  is  suitable  to  dynamic  programming  and  its 
computational  efficiency  is  discussed  in  Section  2.5. 

Section  2.6  is  concerned  with  the  state  dependent  model  where  the 
state  of  the  system  is  defined  to  be  the  number  of  previous  repairs 
performed  on  that  system.  In  other  words,  after  each  repair,  the  new 
failure  rate  will  depend  on  the  number  of  previous  repairs  caused  by 
either  failure  or  planned  preventive  maintenance  operation.  A  review  of 
this  model,  which  was  originally  developed  by  D.G.  Nguyen  and  Murthy 
(1981),  along  with  assumptions,  uniqueness  and  existence  of  the  optimal 
solution  and  a  computational  algorithm  are  described  in  this  section. 
This  computational  algorithm  becomes  inefficient  when  the  number  of 
time  intervals  in  each  renewal  cycle  is  large.  Lemma  2.1  and  policy  I' 
will  permit  us  to  develop  a  heuristic  procedure  which  is  computationally 
more  efficient.  This  is  achieved  by  finding  a  good  estimate  for  N*, 
the  optimum  number  of  planned  repairs  in  each  renewal  cycle.  A 
numerical  example  for  systems  with  Weibull  lifetime  distribution  and  an 
application  to  lemma  2.1  are  also  presented. 

Our  work  extends  the  maintenance  models  to  an  age-dependent  case. 
In  other  words,  the  amount  by  which  recovery  after  a  repair  is  incom- 
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plete  is  related  to  the  length  of  previous  operating  time.  Section  2.7 
contains  two  suggestions  on  how  that  age  factor  can  be  incorporated 
into  the  analysis  of  policy  I. 

1.  Since  the  system's  age  at  the  time  of  repair  is  a  random 
variable  and,  consequently,  it  will  complicate  the  derivation  of  the 
long-term  cost  per  unit  time,  a  simulation  approach  seems  to  be  more 
efficient  to  optimize  the  age-dependent  type  I  policy.  When  age  factor 
after  the  i1^1  repair  is  considered  to  have  the  form 

6.  =  1  +  €{age  at  the  time  of  the  i^1  repair} 

where  is  a  constant  age  deterioration  factor,  and  it  measures  how  fast 

the  system  deterioriates  with  respect  to  its  age. 

2.  To  have  a  model  that  will  be  analytically  tractable,  we  may 
approximate  the  age  by  its  expected  value.  To  be  specific,  the  age 
factor  in  this  case  at  the  time  of  the  i^  preventive  maintenance  action 
can  be  defined  as 

=  1  +  6  {expected  age  at  the  time  of  the  i1^1  repair} 

where  the  age  deterioration  factor,  6,  is  a  constant  parameter  which 
determines  the  relation  between  the  expected  age  and  the  deterioration 
level  of  the  system. 

Section  2.8  deals  only  with  the  simulation  approach.  A  general 
type  I  policy  is  reformulated  such  that  the  resulting  expected  long-term 
average  cost  per  unit  time  is  consistent  with  the  previous  formulation 
derived  in  Section  2.3.  This  mathematical  relation  is  used  to  generate 
the  cost  function  by  means  of  random  numbers.  Among  various  models 
of  this  strategy,  simulation  is  made  of  a  case  in  which  the  age  at  the 
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time  of  repair  will  determine  the  new  failure  rate.  Also,  its  optimization 
with  respect  to  the  expected  long-term  average  cost  criterion  is  dis¬ 
cussed.  Efficiency  and  reduction  in  computation  time  are  obtained  by 
defining  a  relation  among  T1,T2, . •  •  ,Tjj  and  by  estimating  the  optimum 
number  of  preventive  maintenance  operations  in  each  renewal  cycle. 

In  Section  2.9,  we  develop  an  approximate  analytical  model  for  the 
age-dependent  type  I  policy.  The  approximation  is  done  by  replacing 
the  system's  age,  at  the  time  of  the  i ^  repair,  by  its  expected  value. 
A  system  of  nonlinear  equations  which  satisfy  the  necessary  conditions 
of  optimality  will  be  derived .  A  procedure  to  compute  gradients  is 
developed  so  that  the  optimum  solutions  can  be  found  by  applying  a 
numerical  search  technique.  It  is  also  shown  that,  when  underlying 
lifetime  distribution  is  Weibull,  computation  of  gradients  can  be  more 
simplified. 

In  Section  2.10,  policy  I  is  formulated  in  terms  of  steady  state 
availability . 

Chapter  3  is  entirely  devoted  to  policy  II  and  its  optimalities. 

Section  3.2  contains  an  introduction  to  policy  II  by  which  planned 

major  repairs  are  performed  at  T . ;  i  =  1,2, ...  ,N  with  minor  repairs  at 

failures  without  disturbing  the  failure  rates.  The  replacement  takes 
N 

place  at  tM  =  IT.. 

N  i=l  1 

Different  assumptions  regarding  the  hazard  rate  allows  us  to  define 
two  general  models  in  Section  3.3.  "Model  A"  is  applied  when  the 
failure  rate,  after  each  repair,  returns  to  the  same  value  as  at  the  time 
of  replacement  but  with  a  larger  slope  than  if  replacement  had  taken 
place.  "Model  B"  deals  with  systems  for  which,  after  each  repair,  the 
failure  rate  declines  with  unchanged  slope,  but  it  does  not  necessarily 
return  to  the  point  it  started  from. 
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In  Section  3.4,  a  general  mathematical  representation  of  the 

expected  long-term  average  cost  per  unit  time  is  derived  for  policy  II. 

Optimal  policy  II  is  the  solution  to  a  nonlinear  programming  prob¬ 
lem  (NPII)  which  is  found  in  Section  3.5.  First  three  special  cases  are 
discussed . 

1.  Replacement  at  time  T  (constant)  with  minor  repairs  at 
failures.  A  solution  to  this  problem  for  Weibull  lifetime  distributions  is 
also  obtained. 

2.  The  system  ages  rapidly. 

3.  The  system  ages  slowly. 

In  section  3.6,  a  dynamic  programming  formulation  to  type  II 

policy  is  defined  and  its  numerical  computational  efficiency  is  described. 

To  ensure  the  convergence  of  various  optimum  search  techniques 
applied  to  this  problem,  we  investigate  the  existence  of  the  minimum 
point  in  Section  3.7.  Along  with  defining  some  notations,  relations 
satisfying  the  necessary  condition  of  the  optimality  are  also  derived. 

Theorem  3.1  discusses  the  conditions  under  which  the  minimum  point 
can  exist.  Proper  initial  values  for  searching  are  obtained  by  using 
Theorem  3.2. 

Section  3.8  contains  a  review  of  the  state-dependent  model  devel¬ 
oped  by  Nguyen  and  Murthy  (1981).  This  model  assumes  that  the 
number  of  repairs  affects  the  future  failure  rates. 

i 

Section  3.9  defines  two  forms  of  failure  rates  by  which  the  j 

system's  age  can  be  incorporated  into  the  type  II  policy.  One  is  called 

the  A-type  model  with  the  age  factor  represented  as  I 

i 

.  I 

1  j 

9i  =  1  +  IT.  (1.3) 

1  j=l  J 
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The  other  is  the  B-type  model  with  the  age  factor  represented  as 

i 

6.  =  IT. 

1  i=l  } 

Section  3.10  is  devoted  to  the  optimality  of  A-type  age-dependent 
policy  II.  Necessary  conditions  for  optimality  are  derived  for  both 
general  and  Weibull  distributions.  Theorem  3.3  describes  the  order  of 
the  optimal  scheduled  repair  times  when  the  underlying  lifetime  distri¬ 
butions  have  the  Weibull  form.  Other  properties  of  the  optimal  solution 
are  exemplified. 

The  age-dependent  policy  II  of  type  B  and  its  optimality  is  found 

in  Section  3.11.  The  analysis  involves  the  formulation,  necessary  and 

sufficient  conditions  of  the  optimality.  It  is  also  shown  that  for  a  fixed 

value  of  N,  the  optimal  repair  time  is  a  solution  to  an  integral  equation. 

Theorem  3.4  proves  a  useful  inequality  relating  the  optimum  scheduled 

repair  times  and  the  number  of  major  repairs  in  each  renewal  cycle. 

This  theorem  is  then  used  to  find  a  simpler  form  for  the  sufficient 

conditions  when  the  underlying  lifetime  distribution  has  the  Weibull 

* 

form.  The  optimum  number  of  the  planned  repair  times  N  ,  is  computed 

in  terms  of  optimum  scheduled  repair  times,  and  it  is  used  to  obtain  a 

* 

lower  bound  for  N  in  terms  of  maintenance  costs  and  the  failure 
distribution  parameters  (e.g.  Weibull).  Two  analytical  solutions  are 
derived  for  two  cases  of  Weibull  failure  distribution. 

In  Section  3.12  the  policy  II  is  formulated  in  terms  of  steady  state 
availability. 

The  major  results  and  general  conclusions  are  summarized  in 


Chapter  4. 


CHAPTER  2 


POLICY  I 

SCHEDULING  MAXIMUM  TIMES 
BETWEEN  REPAIRS 

2.1  INTRODUCTION 

This  chapter  deals  with  the  optimality  and  the  analysis  of  policy  I 
by  which  a  repair  is  carried  out  at  failure  or  at  a  planned  repair  time, 
whichever  occurs  first.  The  emphasis  is  on  the  optimum  schedules  for 
maintenance  or  repairs,  in  order  to  minimize  long-term  average  operating 
cost  per  unit  time.  This  policy  is  defined  in  Section  2.2.  A  mathe¬ 
matical  function  representing  the  previous  history  of  the  system  and  an 
expression  for  long-term  average  cost  per  unit  time  are  found  in 
Section  2.3.  The  problem  is  formulated  as  nonlinear  programming  and 
dynamic  programming  in  Sections  2.4  and  2.5,  respectively.  Existing 
results  on  state- dependent  models  are  reviewed  and  a  related  improved 
computational  algorithm  is  suggested  in  Section  2.6.  Sections  2. 7-2.9 
are  devoted  to  the  age-dependent  model  and  its  optimality  by  means  of 
a  simulation  technique  and  a  suggested  approximation  procedure.  The 
problem  is  formulated  in  Section  2.10  in  terms  of  steady  state  availability. 

2.2  POLICY  I 

Simple  equipment  such  as  a  car  engine  may  break  down  from  time 
to  time,  needing  repair  to  be  restored  to  operation.  Usually  repair 
actions  for  a  system  having  an  increasing  (or  possibly  nondecreasing) 
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failure  rate  cannot  last  forever,  and  at  some  point  in  time,  replacement 
becomes  inevitable.  It  is  perhaps  justifiable  to  perform  the  repairs 
prior  to  failures.  To  analyze  this  problem,  a  strategy  called  policy  I 
has  received  the  most  attention  in  the  literature.  The  following  is  a 
generalization  of  the  type  I  policy. 

Policy  I.  The  system  is  repaired  for  the  i^1  time  in  the 
case  of  failure  or  after  Tj  time  units  of  continuing  operation 
(from  the  last  maintenance  action),  whichever  occurs  first 
(i<N).  It  is  replaced  when  i=N.  This  process  will  be  re¬ 
peated  indefinitely. 

Using  this  policy  contributes  to  the  reduction  of  irregularity  in 
breakdowns  by  performing  planned  repairs  before  the  system  fails.  A 
discussion  on  this  policy  in  a  practical  case  will  depend  on  considera¬ 
tions  not  all  of  which  can  be  incorporated  into  a  mathematical  analysis. 
Therefore,  the  following  general  assumptions  will  be  made. 

1.  Repairs  are  carried  out  instantaneously  and  there  will  be  no 
queuing  problem. 

2.  Information  about  the  state  of  the  system  (failed,  not  failed)  is 
obtained  without  any  inspection.  In  other  words,  the  system's  state  is 
observable . 

3.  The  failure  time  probability  density  function  (p.d.f.)  is  known 
after  each  repair  action  and  when  the  system  is  new. 

4.  The  system  is  subject  to  stochastic  failures. 

5.  Distributions  of  failure  times  after  each  repair  are  not  neces¬ 
sarily  identical. 

We  will  make  additional  assumptions  as  the  parameters  of  this  model 


are  defined. 
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Several  different  criteria  have  been  considered  to  determine  the 
optimum  planned  repair  times,  among  which  minimizing  the  long-term 
average  cost  and  maximizing  the  limiting  efficiency  (availability)  are 
widely  used.  From  the  optimization  point  of  view,  these  two  criteria 
are  actually  equivalent  (see  Section  2.10). 

We  consider  an  optimum  policy  I  for  which  the  long-term  average 
operating  cost  per  unit  of  time  is  minimized.  Each  cycle,  which  is 
defined  to  be  the  time  between  two  successive  replacements,  consists  of 
N  time  periods  (T\,  i=l,...,N),  i.e.  T.  is  the  time  between  the  (i-l)st 
and  i^1  repair  actions. 

In  this  model  the  maintenance  costs  consist  of  replacement  cost 
CR,  repair  cost  C  ,  and  breakdown  cost  CB  for  each  failure.  These 
costs  are  assumed  to  be  constant,  even  though,  in  practice,  they  might 
well  be  functions  of  age,  number  of  repairs,  and  other  external  factors. 

If  CQ,  CR,  and  especially  CR  are  highly  sensitive  to  the  inflation 
rate  and  other  economic  factors,  then  long-term  total  discounted  cost 
criterion  seems  to  be  preferable.  In  most  situations,  the  costs  can  be 
interpreted  as  expected  values  of  corresponding  random  variables. 

CQ  represents  all  costs  caused  by  performing  a  repair,  including 
the  cost  associated  with  a  constant  repair  time. 

Breakdown  cost  CR  consists  of  all  additional  costs  incurred  due  to 
failure.  When  planned  repair  is  carried  out,  only  CQ  is  incurred. 

It  is  also  possible  to  define  CD,  CD,  and  C  as  functions  of  time 
to  adjust  with  depreciation  or  wear.  But  as  we  said  earlier,  it  is 
assumed  here  that  they  are  constant. 

In  the  next  section,  a  general  type  I  policy  is  formulated. 
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2.3  POLICY  I  FORMULATION 

A  component  or  a  system  subject  to  stochastic  failure  can  be 
described  by  its  failure  (hazard)  rate  r(t) 

r(t)dt  =  Pr[t<x<t+dt|x>t]  (2.1) 

where  x  is  a  random  variable  specifying  the  failure  time.  The  cor¬ 
responding  cumulative  distribution  function  is 

F(t)  =  1  -  exp[-J*  r(s)ds]  (2.2) 

Let  =  gi(T1,T2, . . .  ,T._1;ti,t2, .  • .  ^j.j)  be  a  positive  real 
function  to  describe  the  past  aging  of  the  system  up  to  the  (i-l)st 
preventive  maintenance  action.  In  this  expression,  Tj  (constant 
planned  repair  time)  and  x.  (failure  time)  are  measured  from  the  most 
recent  repair  or  preventive  maintenance.  This  allows  us  to  formulate 
the  type  I  policy  so  that  the  failure  rate  is  a  function  of  age  of  the 
system  and/or  number  of  repairs  performed. 

Failure  characteristics  of  the  system  after  the  (i-l)st  repair  (or 
replacement  if  i=l)  can  now  be  described  by  the  conditional  cumulative 
failure  distribution  function  Fj(t|8._^)  and  the  corresponding  conditional 
hazard  rate  r.(t|8._^),  which  is  assumed  to  be  strictly  increasing  in 
Ti,T2,...,TN,  t,  and  also  in  the  number  of  repairs  (i-1).  Assume 
that  F.(t|8._p  and  r.(t|8._^)  are  both  continuous  and  differentiable  for 
any  positive  value  of  the  vector  T  =  (TWT2, . . .  ,TN)  and  t. 

In  order  to  have  a  meaningful  preventive  maintenance  policy,  we 
assume  that  after  each  repair,  the  system  does  not  fully  recover,  and  it 
is  more  likely  to  break  down  than  if  a  replacement  had  taken  place 
("repaired  not  as  good  as  new").  This  is  achieved  by  making  the 
failure  rate  function  after  any  repair  satisfy 

ri(tlei-1)  *  ri-l(t|ei-2)  *  •*•  *  riWo) 

*— -  — . . m - 


(2.3) 
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Figure  2.1- Various  Linearly  increasing 
Failure  Rates  in  a  cycle. 

The  function  r.  might  depend  on  i-1,  TM,T._2 . Tx.  Figure  2.1 

illustrates  this  fact  when  the  failure  rate  is  linearly  increasing.  Note 
that  repairs  will  reduce  the  failure  rates  considerably  such  that 

ri(°|0i-l)  =  ri-i(°!0i-2)  =  ‘  ••  58  ri(0)  (2.4) 

During  the  ith  period,  the  probability  of  a  failure  before  the  next 
preventive  maintenance,  given  that  (i-1)  repairs  have  been  performed, 
is  expressed  by  Fj(t.  |6j_^),  and  the  expected  breakdown  cost  for  a 
cycle  consisting  of  N  period  is 


N 

Total  Breakdown  Cost  =  Cn  I  F.(T.|6.  ,) 

Bi=l  1  1  1-1 


(2.5) 


The  expected  total  cost  per  cycle  denoted  by  R(N,T)  is  found  to  be 


R(N,T)  * 


CR  +  <N-DC0 


N 


+  C 


B 


(2.6) 
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where 

T  =  (TlfT2...,TN) 

To  calculate  the  expected  length  of  a  cycle,  we  define  the  following: 

tj  =  Random  variable  with  p.d.f.  f(t|Q._^>  denoting  the 
failure  time  in  period  i,  measured  immediately  after  the 

ct 

(i-1)  repair  (or  replacement  for  i=l). 

Tj  =  Scheduled  repair  time  measured  immediately  after  the 

ct 

(i-1)  repair  (or  replacement  for  i=l). 

X.  =  Age  of  period  i,  which  is  either  Ti  or  t^,  whichever 
occurs  first. 

Thus  we  have 


Xj  =  min(ti,T.) 


(2.7) 


l 


using  this  notation ,  one  possible  age  function  is  0.  =1+61  X. ;  where 

1  j=l  J 

>0  is  an  aging  factor.  The  expected  length  of  period  i  is  computed  by 

T. 

y.  $  E(X.)  =  /olsf.(s|0..1)ds  +  TiJ^f.(s|e._1)ds 
Making  a  change  of  variables 


(2.8) 


s  =  U  ;  fi(sl8._1)ds  =  dV 
or 

ds  =  dU  ;  Fi(s|6i_1)  =  V 


and  using  integration  by  parts,  we  get 


T. 


y(  =  -  /0lFi(s|e|_1)ds  *  TjFjCTjie^j) 
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or  yt  =  Tj^Tjie^)  ♦  PjCTjie^)]  -  J^FjCsie^ds 
=  Tj  -  s}  Fi(s|ei_1)ds 

=  4ids  -  -ToiFi(siei-1)ds 

and  finally 

yj  =  F^siej.pds 


\ 

(2.9)  i 


Therefore,  the  expected  length  of  a  cycle  denoted  by  L(N,T)  is  | 

expressed  by 

! 


.  N  NT. 

L(N,T)  ^  E(  2  X.)  =  I  f1  F.  (s 1 0-  i  )ds  (2.10) 

i=l  1  i=l  0  1  11 

According  to  the  renewal  reward  theorem  [Ross  (1970)],  the 
expected  long-run  average  cost  is  the  expected  total  cost  during  a 
cycle  divided  by  the  expected  length  of  that  cycle.  Thus  the  expected 
long-term  average  cost  per  unit  of  time  for  policy  I  is 


C(N,T)  = 

L(N,T) 


Cr  ♦  (N-DC0  +  CB  I^TjIGj.!) 


N  T._ 

z  41f  (tie.  ,)dt 
i=l  0  i  1  1 


(2.11) 


In  the  next  section,  optimization  of  preventive  maintenance  policy  I 
is  given. 

2.4  OPTIMAL  POLICY  I 

Our  goal  is  to  find  a  set  of  positive  time  periods  (T.,  i  =  1,2, ...  ,N) 
so  that  the  expected  long-run  average  cost  per  unit  time  C(N,T)  is 
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minimi2ed.  This  can  be  expressed  by  the  following  nonlinear  program¬ 
ming  problem  (NPI): 


Find  N*  and  T*  =  (T1(N),T2(N),...,TN(N)) 

to  minimize 


C(N.T)  = 

R(N.T) 

L(N,T) 

(2.12) 

subject  to 

Tj  >  0  i 

«  1.2,. ...N 

(2.13) 

NH 

(2.14) 

N  integer 

(2.15) 

First  we  consider  three  special  cases: 


Case  I.  No  Repairs  (N=l) 

This  case  deals  only  with  replacements,  and  repairs  are  not  carried 
out.  If  N=l,  the  C(N,T)  will  reduce  to 


C(l,Ti) 


CR  +  C^tCTQ 
J^F^Odt 


(2.16) 


By  setting  its  first  derivative  equal  to  zero,  we  obtain  the  first 
necessary  condition  for  a  minimum 


Ti_ 

CBfi(Tx);0  FxCOdt  -  Fx(T!)[Cr  +  CgFiCTx)]  =  0  (2.17) 


fj(Ti) 

Dividing  both  sides  of  (2.17)  by  F^T^  and  replacing  -  with 

Fi(Tx) 


fi(Ti)  yields: 
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A  T1 w  CR 

GxCTO  S  r^TJf  F i(t)dt  -  Fx(Tx)  =  ^  (2.18) 

o  uB 

An  optimum  type  I  policy  is  obtained  by  Tx(l),  which  satisfies 
equation  (2.18),  and  the  corresponding  minimum  expected  cost  value  is 

C*(l,Tx(l))  =  CBrx(Tx(l))  (2.19) 

Barlow  and  Proschan  (1965)  and  Barlow  and  Hunter  (1959)  have 
considered  a  similar  model  in  which  they  proved  that  optimal  solution 
Tx(l)  may  exist.  Their  discussion  follows  from  the  fact  that  the  left- 
hand  side  of  equation  (2.18)  is  increasing  in  Tx  provided  hazard  rate 
rx(t)  is  increasing  in  t.  Mathematically: 


dGx(Tx) 

dTx 


T 

r’i(ti)/0  1  F  x(t)dt  +  rl(Tx)F1(T1)  -  fx(Tx) 


=  Fx(Tx) 


.  Tx 

ri(Ti)/0  F  x(t)dt 


Fi(Tx) 


ri(Tx) 


fi(Ti) 

Fi(Tx) 


T:_ 

=  r’x(Tx);0  F  x(t)dt  £  0 


The  right-hand  side  of  equation  (2.18)  is  constant  for  any  value 
of  Tx<  We  conclude  that  C(1,T)  has  at  most  one  minimum  point.  But 
from  (2.15) 


C(l,Tx)  -*■  ®  as  Tx  -*•  0 

and 


C(l.Tx)  - 


CR+CB 

Mi 


as  Tx  -*■  oo 
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where  ^  is  the  expected  failure  time.  It  was  assumed  previously  that 
C(l,Tt)  is  continuous.  The  two  possibilities  are  shown  in  Figures 
2.2(a)  and  2.2(b). 

In  Figure  2.2(a),  clearly,  Tx(l)  =  »,  but  in  2.2(b),  there  exists  a 
finite  positive  Tx(l)  for  which  C(1,T1(1))  <  C(l,»). 


Case  II.  Rapid  Aging 

In  this  case,  the  equipment  deteriorates  very  rapidly  and/or 
repairs  are  destructive.  For  any  age  function  9i  >  1;  (i£l)  and  0Q  =  1, 
the  conditional  failure  rate  may  be  defined  as 


ri(t|ei-l}  =  ei-lr(t)  (2.20) 

from  which  the  related  conditional  cumulative  distribution  function  is 
expressed  by 

ri<tl6i-l)  =  1  '  expt-0i-l4  ri^s)dsl 

t  6i-l 

=  1  -  expl-j*  r.(s)ds]  (2.21) 

If  0j_j  (i  =  2,3,...,N)  is  large,  then 

Lim  F.(t|0.  ,)  =  1  i  =  2,3 . N 

Vl~ 

Clearly  no  planned  repair  must  be  scheduled,  and  the  resulting 
expected  long-term  average  cost  per  unit  time  given  by  (2.11)  becomes 


C(l.T)  = 


CR  +  CrFj(Tj) 


B‘ 

"tt: — 

x0  r i(t)dt 


which  is  the  result  of  Case  I. 
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Case  III.  Slow  aging 

This  situation  deals  with  the  case  in  which,  after  each  repair,  the 
system  becomes  as  good  as  new.  Mathematically,  this  is  described  by 

ri(tiei_1)  =  r(t)  (2.22) 

Combining  this  relation  and  equation  (2.11),  we  get 
CR  ♦  (N-1)C0  *  CB  I F(T|) 

C(N,T)  =  - — -  (2.23) 

N  T. 

*  L  F(t)dt 
i=l  0 

Taking  its  derivatives  and  equating  them  to  2ero  gives 

NT.  N 

CBf(Tj)  F (t)dt  -  F(T.)[Cr  +  (N-1)Cq  +  Cr  I^F(T.)]  =  0 

for  j  =  1,2,. . .  ,N  (2.24) 

which  can  be  simplified  to 

NT.  N  Cc  +  (N-l)C 

r(T.)  I  /  1  F (t)dt  -  I  F(T.)  =  -5— r - -  (2.25) 

]  i=l  °  i=l  1 

The  minimum  expected  cost  per  unit  of  time  is  obtained  by  dividing 
NT. 

(2.24)  by  F(T.)  Z  J  1  F(t)dt 
i  i=l° 

C(N,T  )  =  Cfir(T.(N))  for  all  j  and  K  (2.26) 

which  means 

r(T.(N))  =  r(TK(N))  for  all  j  and  K  (2.27) 

Since  the  failure  rate  r(t)  is  assumed  to  be  an  increasing  function  of  t, 

we  have 

Tj(N)  =  TK(N)  =  T  for  j  and  K  (2.28) 
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And  now  equation  (2.25)  can  be  written  as 


T  CR  +  (N-l)C 

Nr(T)/T(t)dt  -  NF(T)  =  — - -p - -  (2.29) 

o 

Therefore,  the  optimum  repair  time  for  each  period  T  is  the  solution  to 

T  c  +  (N-l)C 

G(T)  £  r(T)J0  F(t)dt  -  F(T)  =  -  -  h(N) 

B 

(2.30) 

The  right-hand  side  of  (2.30)  is  increasing  in  N  if 


+  (N-1)Cq 
NC^ 


CR  +  NCo 
(N+i)cB 


which  implies 


CR  >  Co 


(2.31) 


(2.32) 


Obviously  if  CR  <  CQ,  then  h(N)  is  a  decreasing  function  of  N. 
Since  the  failure  rate  is  increasing,  the  left-hand  side  of  (2.30)  is  also 
increasing  in  T.  To  show  this,  we  compute 


=  r'(T)/J  F (t)dt  +  r(T)F(T)  -  f(T) 

=  r’(T)/T  F(t)dt  +  •  F(T)  -  f(T) 

0  F  (T) 

=  r’(T)/J  F(t)dt  S  0  (2.33) 

If  replacement  cost  CR  is  greater  than  repair  cost  CQ,  then  h(N) 
decreases  as  we  increase  N,  and  consequently  from  (2.33) 


T(N+1)  <  T(N) 


for  all  N 


(2.34) 
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which  implies 

C(N+l,f )  <  C(N,T)  for  all  N 


(2.35) 


[See  Eq.  (2.26)] 

Therefore,  if  the  replacement  cost  is  more  than  repair  cost,  it  is 
optimal  not  to  replace  the  existing  system  as  long  as  it  is  needed.  The 
optimum  repair  time,  T(N),  is  the  solution  to 


r(T )/J  T (t)dt  -  F(T)  =  gS 

B 


(2.36) 


where  is  the  limit  of  h(N)  as  N  approaches  infinity. 

B 

If  the  replacement  cost  CD  is  less  than  the  repair  cost  C  ,  then 

K  O 

IhN)  is  decreasing  in  N  and 


T(N+1)  <  T(N)  for  all  N 


(2.37) 


and  from  (2.26),  we  conclude 

C(N+1,F)  <  C(N,T)  for  all  N  (2.38) 

Thus,  the  optimum  number  of  scheduled  repairs  in  this  case  is 
zero  and  only  replacements  are  carried  out.  The  procedure  is  again 
the  same  as  for  Case  I  (N=l). 

Note  that  the  breakdown  cost  Cfi  affects  the  above  strategies  by 
changing  the  length  of  T(N)  and,  as  a  result,  the  value  of  C(N,T), 
but  it  does  not  affect  the  optimal  value  of  N.  This  is  because  a  reduc¬ 
tion  in  Cg  will  result  in  a  larger  T(N).  [h(N)  is  decreasing  in  Cg.] 

When  CR  =  CQ,  at  the  time  of  failure  or  planned  repair,  we  can 
either  replace  the  system  or  perform  a  repair. 

Case  III,  which  is  also  a  special  case  of  a  model  developed  by 
Nguyen  and  Murthy  (1981),  does  not  seem  to  be  appropriate  in  prac- 


tice,  because  the  system's  age  and  number  of  repairs  are  ignored.  In 
order  to  have  a  nontrivial  solution,  we  admit  that  the  generalized 
policy  I  is  applied  when  replacement  cost  is  greater  than  repair  cost 

Co- 

Now  we  turn  our  attention  to  the  general  form  of  (NPI)  and  employ 
Dynamic  Programming  to  reformulate  it. 


2.5  DYNAMIC  PROGRAMMING  FORMULATION 

The  objective  is  to  select  a  set  of  positive  repair  times  T  = 
(T1(N),T2(N),...,T3(N))  so  that  the  expected  average  cost  per  unit 
time  given  by 


C(N,T) 


R(NiT) 

L(N.T) 


(2.11) 


is  minimized  where  R  and  L  have  been  defined  by  (2.6)  and  (2.10).  It 
is  well  known  [Barlow  and  Proschan  (1965)]  that  T,  which  minimizes 
C(N,T)  will  also  minimize  a  related  problem  given  by 


-Ti(P)  =  R(N,T)  -  pL(N,T )  (2.39) 

*  _ 

Where  p  =  p,  the  minimum  of  C(N,T)  is  such  that 

JtQ  (P)  =  0  (2.40) 

Concentrating  upon  a  single  cycle,  we  define  the  following  recur¬ 
sive  equation 

-^K(P,9K-i)  =  Min{CK  +-£K+1(MK)} 
aK 


; 


(2.41) 
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CK  is  the  cost  incurred  during  each  period,  and  it  is  equal  to 

T 

CK  =  a-nK)cR*nKc0^cBrK(TKl6K.1)-pj'0KrK(t|9K.1)dt 

for  K  =  1,2 . N  (2.42) 

when 

nK  =  1  if  1  S  K  <  N 

nK  =  0  if  K  =  N 

We  now  define  eK  as  being  proportional  to  the  system's  age  at  the 
time  of  the  most  recent  repair.  This  can  be  written  as  a  transition 
equation 


0K  =  0K-1  +  Min(xK,TK)  K  f  0 


(2.43) 


eo  =  1 


where  €  is  a  multiplier  to  determine  the  deterioration  rate  with  respect 
to  age.  0^,  which  is  a  random  variable,  has  the  following  distribution 


0K  =  6K-1  +€TK 

with  probability 

F(tkI9k.1) 

=  0K_1  +  €tr 

with  probability 

e'<tkI|W 

(2.44) 

The  boundary  conditions 

to  (2.41)  are 

-^N+1(P  ,0N)  =  0 

and 

(2.45) 

-A(Mn)  =  0 

(2.46) 

In  most  cases  it  seems  unlikely  to  find  an  analytical  solution  for 
this  problem.  This  is  due  to  complexity  caused  by  transition  equation 
(2.43)  and  p.  As  an  example,  consider  a  system  in  which  the  failure 


characteristic  follows  the  Weibull  distribution.  In  this  situation,  finding 
the  distribution  of  0^  requires  numerical  calculations,  p*  is  found  at  the 
last  stage  and  its  different  values  must  be  kept  at  each  stage.  It  is 
most  likely  that  carrying  out  these  calculations  will  require  a  lot  of 
execution  time  and  storage  space,  and  the  results  will  lack  precision. 

2.6  STATE  DEPENDENT  MODEL 

We  define  the  state  of  the  system  to  be  the  number  of  repairs 
performed.  A  mathematical  representation  of  this  model  is 

ri(t|0i_i)  =  r.(t)  (2.46) 

where  i  represents  the  number  of  repairs .  This  case  describes  a 
situation  when  the  failure  rate  function  increases  with  the  number  of 
repairs  as  well  as  time  since  the  most  recent  repair,  but  the  system's 
age  at  the  time  of  repair  does  not  play  any  significant  role  in  its  future 
performance.  This  problem  has  been  studied  by  D.G.  Nguyen  and 
D.N.P.  Murthy  (1981). 

From  equation  (2.11),  the  expected  cost  per  unit  time  for  this  case 
is 

CR  ♦  (N-1)C0  -  CB  I  F.(T.) 

C(N,T)  =  - rT: - — -  (2.47) 

Z  /  1  F.(t)dt 
i=l  0  1 

As  Nguyen  and  Murthy  (1981)  have  shown  by  differentiating 
(2.47)  with  respect  to  the  T.  (i  =  1,2, . . .  ,N),  the  optimal  preventive 
maintenance  ages  satisfy 
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They  have  proved  that  under  the  following  assumptions 

i)  r.(t)  is  strictly  increasng  to  infinity  (in  t) 

ii)  r.+1(t)  *  r4(t),  t  >  0 

iii)  r.+1(0)  =  r.(0) 
the  following  are  true: 

1.  For  a  fixed  value  of  N,  the  optimal  policy  I  exists  and  T^N); 

(i  =  1,2,...,N)  is  finite,  unique  and  decreasing  in  i. 

* 

2.  There  exists  an  optimal  N  for  which 

min[N|AC(N,f  2  0]  £  N*  £  max[N|AC(N,T)  S  0] 

N  N 


where  AC(N,T)  =  C(N+1,T)  -  C(N,T) 


Nguyen  and  Murthy  also  proposed  a  computation  algorithm  as 
follows : 

i)  Set  N  =  1 

ii)  Obtain  TjCN);  i  =  using  equations  (2.48)  and 

(2.49) 

iii)  If  T(N)  2  T(N-l),  got  to  step  (v) 

iv)  Set  N  =  N+l  and  go  to  step  (ii) 

*  *  _ 

v)  N  =  N-l;  compute  C(N  ,T)  from  equation  (2.47) 

(Note:  step  (iii)  is  omitted  for  N  =  1.) 

They  have  suggested  that  step  (ii)  be  carried  out  by  solving 
identities  (2.48)  and  (2.49)  simultaneously.  However,  applying  any 
optimum  seeking  technique,  such  as  gradient  method,  directly  to  equation 
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(2.47)  also  provides  satisfactory  results.  Either  technique  seems  to 

require  numerical  integration  when  underlying  lifetime  distributions  are 

assumed  to  have  the  Weibull  form.  Consequently  for  a  large  N,  the 

* 

above  algorithm  is  not  likely  to  be  efficient.  For  instance,  if  N  =49, 

then  the  last  iteration  of  the  algorithm  alone  requires  (50)(.2)(15)/60 

=  2.5  minutes  of  computer  time  for  numerical  integration  (assuming  that 

each  numerical  integration  needs  .2  seconds  and  step  (ii)  will  converge 

after  15  iterations).  In  other  words,  to  carry  out  this  algorithm  when 
* 

N  =49,  approximately  one  hour  of  computer  time  is  spent  to  integrate 
numerically. 

* 

The  following  lemma  will  help  us  find  an  approximate  value  for  N  , 
which  will  be  useful  in  later  numerical  calculations. 


Lemma  2.1 

For  constant  (i  =  1,2 . n),  a  and  b  (b  >  0),  the  integer 

function 


P 


n 


a  +  bn 
n 


n  =  1,2,...  (2.50) 


has  at  least  one  global  minimum  if 

Mi  *  P2  *  •••  *  Mn  i  •••  >  0  (2.51) 

Proof.  Let  Apn  -  pn+1  -  pn  and  from  relation  (2.50)  we  have 

n 

-%i«  *  b<  I/i  -  nVi> 

Ap  =  - 3 - — - 

n+1  n 

(  I  Mj)(  I  M;) 
i=l  1  i=l  1 


(2.52) 
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Note  that  if  a  >  0,  and  p.+1  -  pi  for  all  i,  then  Apn  <  0,  which 

# 

implies  that  p  decreases  as  n  increases  or  n  =  ».  If  a  <  0,  then 
n  n 

Apn  >  0  for  all  n.  Since  I  p^  -  npn+^  >  0  and  p.  £  pn+1;  i  =  l#2,...,n 

(from  2.52)  and  consequently  n  =1.  For  a  >  0  and  b  >  0,  first  we 
will  show  that  if  Apn  >  0,  then  Apft+1  >  0.  Identity  (2.52)  is  positive 
when 

-apn+i  ♦  b(^p.  ‘  n'Jn+l)  >  0  (2‘53) 


or 


n 

.2  Mi 
i=l  1 


n>Jn+l 


^n+1 


a 

B 


Similarly  ApN+1  >  0  if  and  only  if 


n+1 


I  M,  -  Cntl)Mn+2 

-  > 


^n+2 


which  can  also  be  written  as 


(2.54) 


(2.55) 


n 

I  P;  -  np 
i=l  1 


n+2  +  Mn+1 


^n+2 


n+2 


a 

B 


(2.56) 


We  can  increase  the  left-hand  side  of  inequality  (2.54)  via  re¬ 
placing  pn+1  by  pn+2  and  adding  pR+1  -  pn+2  >  0  to  the  numerator. 
This  will  result  in  inequality  (2.56)  and,  hence  we  conclude  that  Apn  >  0 
implies  0pn+^  >  0. 

Now  we  have  to  prove  that  Apn+j  <  0  implies  Apn  <  0.  For  Apn+1  <  0, 


it  can  be  shown  that 
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n 


I 

i=l 


n^n+2  +  >Vl  '  Mn+2 

*^n+2 


< 


a 

E 


(2.57) 


The  inequality  holds  if  we  decrease  the  left-hand  side.  Replacing 


Mn+2  by  Mn+1  yields 


£  Mi  *  niVl 


a 

E 


(2.58) 


which  means  Apn  <  0  and  the  proof  is  completed. 

One  may  define  the  following  policy: 

Policy  I1:  The  system  is  repaired  for  the  i^  time  in  the  case 

of  failure  if  i  <  N'.  It  is  replaced  when  i  =  N'. 

In  fact  this  strategy  is  the  same  as  policy  I,  but  planned  repair 
times  have  been  selected  so  large  that  failures  occur  prior  to  the 
scheduled  preventive  maintenance.  The  long-term  average  cost  per  unit 
time  is  obtained  by  letting  T-  -*■  ®  (i  =  1,2,...,N)  in  identity  (2.11) 


C(N\») 


CR  “ 


Co  +  <Co+CB>N’ 


N' 

2  M 
i=l 


l 


(2.59) 


where  is  the  expected  failure  time  in  the  i111  period.  It  has  been 
shown  by  Nguyen  and  Murthy  that  for  a  state-dependent  model  with  the 
assumed  failure  rate  function  properties,  p^  §  p.+1  >0;i  =  l,2,....  In 
fact,  this  attribute  exists  in  age  and  state-age-dependent  models.  This 
follows  from  the  fact  that  whenever  r^tle.  ^)  is  less  than  (t 1 0.),  one 
can  see  easily  that  F.(t|0. _^)  <  F^tie.  ^)  and  this  leads  to  Mj+1  <  Mj- 
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Setting  CR-CQ  =  a,  CQ+CB  =  b,  and  N*  =  n,  we  get  equation 
(2.50)  given  in  lemma  2.1,  and  hence  we  conclude  that  optimal  policy  I* 
exists.  There  is  a  unique  N'  which  minimizes  the  expected  cost  per 
unit  time  (2.59)  except  the  case  in  which  two  or  more  values  of  N'  will 
result  in  the  same  amount  of  cost  C(N',»).  In  other  words,  function 
(2.59)  is  either  convex  or  strictly  convex  function  of  N'.  It  is  important 
to  note  that  the  behavior  of  C(N,»),  when  =  p-  (policy  I:),  and 
special  case  III  (policy  I)  are  exactly  the  same.  In  both  models, 
N  =  N'  =od  when  CD  >  C  ,  and  for  CD  <  <2,  no  repair  should  be 

K  O  K  O 

scheduled  (N*  =  N'*  =  1). 

Calculation  of  N1  is  simple.  T.'s  have  been  eliminated,  and  there 

is  no  numerical  integration  to  evaluate.  In  most  cases,  it  can  be  done 

*  * 
by  hand.  But  the  question  is  how  well  N  is  approximated  by  N'  . 

Our  simulation  study  shows  that  N'  is  a  good  estimate  for  N,  and 

due  to  the  lack  of  preventive  maintenance  actions,  the  average  cost  for 

*  * 

strategy  I'  is  higher  than  policy  I;  consequently,  N'  £  N  .  This 
* 

allows  us  to  use  N'  as  a  starting  value  for  N.  Employing  this  heuristic 
procedure  will  save  a  large  amount  of  computer  time  by  avoiding  the 
computation  used  to  evaluate  the  iterations  1,2, . . .  ,N'-1. 

Numerical  Example  2.1 

Referring  to  the  numerical  example  given  by  Nguyen  and  Murthy 
where  CQ  =  5,  CR  =  15,  Cg  =  15,  and  failure  times  follow  the  Weibull 
distributions,  i.e. 

\  a+l 

fj(t)  =  Ajt"  •  exp(-  -i?]-)  i  *  1.2 . N 


f 


1 
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with  constant  shape  parameter  a  =  1  and  scale  parameters  Xi  =  1, 

=  (1.5)1"1;  i  =  2,3, The  optimum  type  I  policy  is  obtained  by 

Tx(3)  =  0.936,  T2(3)  =  0.624,  T3(3)  =  0.416,  N*  =  3  and  C(N,T*)  = 

* 

28.08.  But  our  heuristic  procedure,  policy  I’,  yields  N'  =  2.  By 
selecting  initially  N  =  2,  we  could  have  avoided  more  complicated  compu¬ 
tations  required  for  N  =  1.  The  results  are  shown  in  Table  2.1. 

_JV _ C(N',«) 

1  33.85 

2  31.06 

3  31.80 

Table  2.1.  Optimal  Policy  I’  (N'*  =  2) 

Numerical  Example  2.2 

This  example  illustrates  a  case  for  which  calculations  of  N'  can 
further  be  simplified.  Suppose  each  failure  (and  consequently  repair) 
affects  the  expected  lifetime  of  the  following  period  such  that 

P,  =  (S^ii)u0  =  (1-  ~)M0  i  =  1,2 . n  (2.60) 

i  n2  °  n2  0 


where  Mi  =  P0-  Note  that  p.  is  decreasing  in  i  and  the  expected  length 
of  the  i^1  time  interval  is  less  than  the  expected  (i-l)st  lifetime. 

The  equation  (2.50),  for  this  example,  becomes 


a  +  bn 


n  n 

2  <*-  ~>^o 


i=l 


M. 

n2 


a  +  bn 
,  n-1  N 

Mo(n'  5rT  ) 


We  may  treat  n  as  a  continuous  variable  and  differentiate  pn  with 
respect  to  n.  Setting  that  derivative  to  zero  yields 

(2a  +  b)n2  -  2bn  -  a  =  0  (2.61) 
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The  positive  solution  to  this  quadratic  equation  is 


b  =  Jb2+ab+2a2 
2a+b 


(2.62) 


Setting  a  =  CR-CQ  and  b  =  (C0+Cg),  we  get 


VCB  *  ^CR<2V3Co>+VCR+Co+CB> 

2CR"Co+CB 


(2.63) 


The  optimum  N'  is  either  [n]  or  [n]+l,  depending  on  the  value  of 
C(N',<»).  [n]  denotes  the  largest  integer  less  than  or  equal  to  n. 


2.7  AGE-DEPENDENT  MODEL 

We  are  now  considering  a  model  that  is  applied  to  situations  in 

which  repairs  do  not  reflect  any  increase  in  the  system's  failure  rate. 

But,  due  to  the  length  of  operating  time,  recovery  is  not  fully  achieved 
by  performing  a  repair.  According  to  our  notation,  becomes  an  age 
factor  that  relates  the  system's  age  at  the  end  of  the  i^  period  to  its 
future  performance.  For  simplicity,  we  assume  that  the  conditional 
hazard  rate  for  each  period  is  proportional  to  the  age  factor  e.^  at  the 
beginning  of  that  period 

ri(t|0i-1)  =  r(t)ei-i  (2.64) 

This  implies  that 

0i  1 

FCt|0i-l)  =  tF(t)]  (2.65) 

and  0._!  is  restricted  to 


e._1  >  1  ,  if  i  >  1 

=  1  ,  if  i  =  1 


(2.66) 
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In  fact,  6.  is  a  random  variable  defined  by 

0J  =  . X.^)  (2.67) 

where  g-  is  an  increasing  function  of  Xx,  X2,...,X._^  (X.  is  given  by 
equation  2.7).  An  appropriate  example  of  g.  is 

i 

0.  =  1+  6 1  X;  (2.68) 

1  j=l  J 

which  satisfies  the  conditions  given  by  (2.66)  and  2.67).  6  is  a 

positive  parameter  which  determines  how  fast  the  system  deteriorates 
with  respect  to  its  age.  A  simulation  study  seems  to  be  appropriate  to 
analyze  this  age-dependent  type  I  policy.  We  may  also  approximate 
relation  (2.68)  by 

i 

0:  =  1+  €ly:  (2.69) 

1  j=l  J 

where 

yj  ^  E(X.)  (2.70) 

This  way,  the  failure  rates  will  depend  on  the  expected  value  of  age  as 

a  constant  rather  than  age  as  a  random  variable. 

In  the  next  two  sections,  the  age-dependent  policy  I  is  formulated 
in  two  ways,  one  via  a  simulation  approach  and  the  other  by  approxi¬ 
mation,  using  equations  (2.69)  and  (2.70)  instead  of  (2.68). 

2.8  SIMULATION  APPROACH 

A  key  step  in  the  analysis  of  the  age-dependent  type  I  policy  via 
simulation  is  to  find  expressions  for  total  cost  per  cycle,  R(N,T),  and 


length  of  a  cycle  L(N,T),  in  terms  of  random  variables  involved.  For 
this  to  be  the  case,  let 


1,  if  failure  time  <  T. 
0,  otherwise 


This  is  equivalent  to 


i  =  1,2 . N 


(2.71) 


Pr{6j=l}  =  F-CT.ISj^) 


(2.72) 


Pr{6.=0}  =  1-Fi(Ti|ei_i) 


(2.73) 


According  to  policy  I,  breakdown  cost  Cg  is  incurred  if  the  system 
fails  to  operate  before  the  scheduled  repair  time.  Given  this  notion, 
and  using  the  definition  given  for  (2.71),  the  total  cost  per  cycle, 
R(N,T)  (random  variable),  is  expressed  by 


R(N.T)  =  CR  +  (N-1)C0  +  CB 


(2.74) 


A  random  variable  representing  the  length  of  a  cycle  is 


L(N,T) 


i=l  1 


=  Z  min(x.,T.) 
i=l  1  1 


(2.75) 


or  by  employing  6.,  it  can  be  shown  as 


L(N,T)  =  M  6.x.  +  (1-6.)T.] 


(2.76) 


-39- 


Let  R^(N,T)  and  L^(N,T)  be  the  result  of  the  i^1  simulation 
experiment  and  let  us  assume  that  n  such  experiments  are  to  be  per¬ 
formed,  where  n  is  a  fairly  large  number  (  50  5  n  £  100).  Then  the 
long-term  average  cost  per  unit  time  consistent  with  relation  (2.11)  is 
obtained  by 

2  R(l)(N,T ) 

C(N,T)  =  — - - -  (2.77) 

I  L(l)(N,T) 


As  one  might  expect,  the  above  formulation  is  a  general  repre¬ 
sentation  of  type  I  policy,  and  it  is  applied  to  state-,  age-,  or  state- 
age-dependent  models.  But  generating  a  random  variable,  x-,  differs 
for  each  case. 


Among  the  most  often  used  random  number  generators  is  the 
multiplicative,  congruential  method.  But  as  Marsaglia  (1968)  has  pointed 
out,  this  technique  provides  unsatisfactory  results  when  sequences  of 
random  vectors  are  to  be  generated.  If  K-tuples  (RN1,RN2, . . .  ,RN^), 
(RH2,RN3, . .  .RNK+1), . . .  of  random  numbers  generated  by  this  method 
are  viewed  as  points  of  a  unit  cube  (K  dimensions),  all  the  points  will 
lie  in  a  relatively  small  number  of  hyperplanes.  This  property  indicates 
that  multiplicative  congruential  method  is  inadequate  for  our  type  I 
policy.  Marsaglia  and  Bray  (1968)  have  suggested  the  use  of  a  com¬ 
posite  generator,  which  mixes  three  congruential  generators:  one  to  fill 
nx  locations,  one  to  select  a  location  from  the  n1(  and  the  third  is  used 
for  good  measure. 


The  following  algorithm  can  be  used  to  generate  n  random  variables 

6.  , 

having  the  distribution  F(tl6._1)  =  l-[F(t)]  ,  which  is  appropriate 

for  the  age-dependent  type  I  policy. 


Step  0 
Step  1 
Step  1 
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Set  0Q  =  1;  i  =  0 
Set  i  =  i+1 

-  6i-l 

Generate  t.,  having  distribution  F(t|0^)  =  l-[F(t)] 
using  random  number  RN. 

Step  3:  If  i  =  N,  stop;  else  go  to  Step  4 
Step  4:  Compute  X.  =  min(t.,T.) 

Step  5:  Compute  0.  =  0^  +  £X. 

Step  6:  Go  to  Step  1 

As  an  example,  consider  a  model  in  which  the  underlying  distri 
bution  is  Weibull,  i.e. 


i-1 

0i  i  =  1  +  I  min(t.  ,T-) 
1-1  j=l  1  1 


i  =  1,2, ... ,N  (2.78) 


(2.79) 


Then  6.  =  1,  otherwise  6.  =  0. 

Another  problem  that  can  play  a  significant  role  in  finding  the 
optimal  policy  I  is  selecting  Tx,  T2,...,TN.  Obviously,  simulating  all 
possible  values  will  lead  to  inefficiency.  However,  it  is  reasonable  to 
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say  that  the  maximum  length  of  period  i  (a  good  system)  should  be 
longer  than  the  one  of  period  (i+1)  (not  as  good  a  system).  Our 
simulation  study  also  proves  this  fact  when  the  lifetime  distribution  is 
Weibull.  Given  this  notion,  it  is  reasonable  to  reduce  the  dimensionality 
of  the  search  by  assuming  a  proportional  relationship  among  the  T^.  We 
may  define  T.+1  =  aT.;  (i  =  1,2,...,N-1;  0<a<l),  and  then  find  the 
optimal  Tlf  N,  and  a.  This  method  is  quite  simple,  but  for  a  large  N 
and  a«l,  there  is  a  K  such  that  for  i>K,  the  value  of  Tj  is  no  longer 
significant.  However,  it  can  also  be  shown  by  numerical  examples  (e.g. 
Weibull)  that  there  is  not  a  constant  proportionality  among  optimal  times 
Tx,T2, . . .  ,TN.  In  the  Weibull  case,  we  obtained  better  results  by 
letting  the  reliability  of  the  system  in  period  i+1  be  related  to  the 
reliability  of  period  i  by  a  constant  value  B,  i.e. 


Pr{t  >  T.+1}  =  [Pr{t  >  T.}] 


B 


Dr  /■foitlWt'9i> 
or  e  dt 


-BJcVlVi) 

e  dt 


(2.82) 


(2.83) 


This  becomes 


Vi(t|8i)dt  ■  B4‘ri(t|8i-i>dt 


T, 


(2.84) 


For  a  Weibull  distribution  with  hazard  rate  r.(t|0._^)  =  A.  e^jt®,  identity 
(2.84)  leads  to 


V°  e.T*?  =  S?  e,  1t?+1 


0+1  i  i+1  "  o+r  i-1  i 


(2.85) 
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or 


1  1 


1 

where  A  =  Ba+^  .  In  most  cases,  the  optimal  A  seems  to  be  close  to  one 
(.95  3  A  £  1). 

Numerical  Example  2.3 

Assume  failure  times  of  a  system  follow  the  Weibull  distribution 

given  by  (2.78)  and  (2.79)  with  \  =  1,  a  =  1,  and  €  =  .2.  The 

maintenance  costs  are  CR  =  15,  CQ  =  5,  and  CR  =  12.  The  optimal 

solution  obtained  by  employing  the  simulation  approach  (by  assuming 

£  £ 

that  the  relation  2.86  satisfying  among  the  Tj)  is  N  =  7,  A  =1, 
T*  =  (0.9  0.83  0.78  0.73  0.71  0.68  0.67)  and  C(7,T*)  =  15.08 
where  C(7,«)  =  18.23.  Figure  2.3  shows  the  behavior  of  C(N,T*)  as  N 
is  incremented.  The  result  is  a  convex  function,  and  it  has  a  global 
minimum . 

In  a  manner  similar  to  the  state  dependent  model  in  Section  2.5, 
the  computations  may  be  reduced  by  finding  an  estimate  for  N.  This  is 
achieved  by  simulating  a  similar  but  not  as  good  policy,  namely  letting 
T  -*•  00  (6.  =  1) -  Equation  (2.74)  becomes 

R(N,T)  =  CR  +  (N-1)Cq  +  NCb  (2.87) 

and  equation  (2.75)  will  change  to 

N 

L(N,f)  =  I  x. 

i  —  1  A 


(2.88) 


because  of  the  higher  cost  in  this  case.  The  optimal  N  (say  N'  )  is 
less  than  the  one  of  actual  policy  I,  and  N'  may  be  used  as  a  starting 
value  for  N.  Figure  2.4  illustrates  the  behavior  of  C(N,»)  as  N  is 
incremented  (for  the  numerical  example  2.3).  It  has  also  a  global 

minimum.  The  estimated  number  of  scheduled  repairs  for  this  case  is 

* 

N'  =4.  By  starting  from  N  =  4,  computations  required  for  N  =  1,  2, 
and  3  could  have  been  replaced  by  fewer  computations  needed  to  obtain 


2.9  APPROXIMATION 

Although  simulation  provides  a  flexible  and  significant  approach  to 
the  age-dependent  type  I  policy,  nevertheless,  one  may  wish  to  have 
rather  an  approximate  analytical  formulation  to  this  problem.  For  this 
to  be  the  case,  we  let 


6.  =  1  +  l  E[min(x.,T.)]  =  1  +  6  Z  y. 

j=l  J  J  j=l  J 


(2.89) 


where  y .  is  the  expected  age  of  the  period  defined  by 


Vi  -  -To1  Wiei.pdt 


(2.90) 


In  other  words,  the  failure  rate  at  the  beginning  of  each  stage 
depends  on  the  expected  current  age  rather  than  the  actual  age. 
Therefore,  the  expected  cost  per  unit  time  (2.11)  can  be  written  as 

N 

CD  +  (N-l)C  +  CD  Z  F(T.  |6.  . ) 
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3yi  - 

Since  =  F (Ti 1 9i_1) /  (2.91)  can  be  rewritten  as 


C(N,f )  = 


N  3y. 

CR  +  (N-1)C0  +  CB  ) 


N 

,iy‘ 


(2.92) 


According  to  our  assumption  r(t |0._^)  =  r(t)6j_1  and  hence 


^f(tl0i-1)  =  exPt_i'olr(t)6i-ldt^ 


T-  0 

=  [exp(-/Qi  r(t)dt)]  i_1 


=  [F  (t)] 


0. 


i-1 


(2.33) 


or 


F(t|0 


F  (t) 


i-1 

1+61  y, 
j=l  ] 


(2.94) 


which  is  an  indication  of  reduction  in  reliability  as  the  expected  age 
increases  and  6  will  determine  the  magnitude  of  that  reduction. 

To  find  an  expression  satisfying  the  necessary  condition  of  the 
optimality,  denote  the  first  partial  derivatives  of  relation  (2.92)  by 


C.(N,f)  = 


and  setting  Cj(N,T)  to  zero  yields 
N  32y.  N  N  3y.  r 

‘CBiIi  ^  i=Vi  '  iWC 


j  =  1.2 . N 


(2.95) 


N 


3T.  i 


R  +  +  CBif1  3TT  > 


=  0 


(2.96) 
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which  can  be  simplified  to 

N  d2yi  N  3y.  N  N  8y;  CR+(N-1)C 

-  Z  /  2  aV  2  Yi  ~  2  (1»t  )  =  - 2  (2.97) 

i=j  9Tj3Ti  i=j  3Tj  i=l  1  i=l  aTi  CB 

j  =  1,2 . N 

Identity  (2.96)  can  also  be  written  as 


Cr+(N-1)C0+Cb 


N  ay. 

2  d-  pr ) 
i=i  aTi 


N 

Z 

i=l 


N 


-C 


a2y. 


B  aT.aTj 


n  dy. 

/  z  ssr  (2.98) 

i-J  3Ti 


a  *  i.2 . n> 


The  left-hand  side  is  just  the  minimum  expected  cost  per  unit 
time.  When  j=N  from  equation  (2.90),  we  have 


3y 


N 


8T.r  F(TN,0N-1) 


kN 


and 


a2y 


N  _ 


aT?T  "  "f(TN,eN-i) 


N 


hence  the  right-hand  side  of  (2.98)  is 


(2.99) 


(2.100) 


RHS  =  -f - SLJLL.  =  cRr(Tw|0N  .)  (2.101) 


Therefore,  the  minimum  expected  cost  per  unit  time  can  be  expressed  by 
C(N,T*)  =  CBr(Tj^|0^_1) 


(2.102) 
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i 

j 

The  above  computations  can  also  be  used  to  simplify  one  of  the 
equations  (j=N)  given  by  (2.97)  to  get 


r(TN!8 


N 

’N-lVi 


N 

-  1  F(T 
i=l 


N|0N-1)  " 


CR  +  (N-l)Cc 


(2.103) 


Evaluating  (2.97  for  j=K  and  K+l  and  subtracting  the  results,  we  get 
another  system  of  equations  for  the  first  necessary  condition  of  the 
optimality. 


32Yi 


K+l3Ti 


N 

/  I 
i=l 


/?  a2yi 

3TK3Ti 


N 

/  1 
i=K 


K  =  1,2, .. . ,N-1 


(2.104) 


Despite  the  fact  that  equations  (2.97)  or  (2.104)  along  with  (2.103) 
provide  a  way  to  obtain  the  optimal  solution,  they  seem  to  be  quite 
tedious  and  difficult  to  solve.  One  might  prefer  to  apply  an  optimum 
search  technique  (e.g.,  conjugate  gradient)  directly  to  identity  (2.91) 
and  have  a  more  convenient  and  efficient  method  for  the  evaluation  of 
optimal  age-dependent  policy  I.  For  this  to  be  the  case,  we  must  find 
expressions  for  gradients  C.(N,T),  j  =  l,2,...,N.  This  in  return 

requires  the  derivations  of 

Rj(N,T)  =  j  =  1,2 . N  (2.105) 

and 

Lj(N ,T)  =  — j  =  1,2, ...  ,N  (2.106) 

To  find  R.(N,T),  it  is  simplest  to  proceed  by  induction.  This  will 


result  in 
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Rj(N.T) 


=  C 


B 


d+  eVy.wTpjrcT.) 

i=l  1  J  L  1 


j-1 

1€  1  yi 
i=l  1 


N  m-1  3y. 


-61(1 


3T. 


m=j+l  i=j 


)|F(Tm) 


m-1 

1 1+  Vl  Y“ 


InFCTJ 


N 

(  I  (  )  =  0  when  j  =  N) 

m=j+l 


j  =  1,2 . N  (2.107) 


y^  is  the  first  term  in  L(N,T)  which  involves  TV.  Hence 


N  dy. 

L.(N,T)  =  1  gTJT 

J  i=j  j 


j  =  1,2 . N  (2.108) 


The  partial  definative  of  yK  with  respect  to  T.  is  obtained  from 
relation  (2.90) 


ayK 

3T 


K-l  3y.  T 


^  =  £(  2  -1)/, 


K 


i  j=i  3Tj 


F  (t) 


K-l 

l+€2  y, 

j=i  J,_- 


lnF  (t)dt  for  i  =  1,2 . K-l 


=  [F 


F<TPj 


i-1 

ll+  €  2  y. 
j=l  J 


=  0 


for  i  =  K 


for  i  =  K,K+1, . . .  ,N 


K  =  1,2,. ..,N 


(2.109) 


In  order  to  represent  the  above  formulas  in  a  rather  simpler 
fashion,  let 


H  -  3yK 
dKi  ~  WT~ 


i  =  1,2, ...N 


(2.110) 
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and  also  let 

K-l 

T  r  -jl+  6  2  y| 

VK  =  -/0  K[^(t)  J  j=1  lnF (t)dt  (2.111) 

K  =  1,2, ... ,N 

Since  lnF(t)  £  0,  then  will  be  positive.  Thus  (2.109)  can  be  re¬ 
written  as 

K-l 

dKi  =  -  £(  2  d..)VK  for  i  =  1,2,..., K-l 

i-1 

=  F(T.)  j=1  for  i  =  K 

•  * 

=  0  for  i  =  K+l , K=2 . N  (2.112) 

As  one  might  expect,  we  can  now  form  a  matrix  d,  whose  elements 
dKi  are  related  by  relation  (2.112).  This  matrix  is  used  for  calculations 
of  both  Ri(N,T)  and  L^N.T). 

Matrix  d  is  shown  by 


11 

0 

0  .  .  . 

0 

‘21 

d22 

0  .  .  . 

0 

31 

• 

d32 

d33.  .  . 

• 

• 

0 

(2.113) 

• 

N1 

dN2 

• 

dN3'  *  ’ 

dNN 

Equation  (2.108)  can  be  calculated  by  using  matrix  d 
N 

L.(N,T)  =  I  d.. 
j=i 


sum  of  column  i 


(2.114) 
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Calculation  of  R.(N,f )  also  requires  partial  summation  along  each 
column  of  matrix  d.  After  obtaining  R.(N,T)  and  L.(N,T),  the 
gradients  are  simply  expressed  by 


C^N.T) 


R.(N,T  )L(N,T)  -  L.(N,T)R(N,T) 
[L(N,T]2 


(2.115) 


As  an  example,  let  us  assume  that  the  underlying  lifetime  distri¬ 
bution  is  Weibull  given  by 


fi(t,ei-l>  =  ^i-l1  exp( - a+i — ) 


(2.116) 


a  >  1  ;  X0,t  >  0 


where 


i-1 

Vi  =  1+  6  1  Y\ 
11  =1  1 


(2.117) 


From  identity  (2.107),  Rj(N,T)  for  this  case  is 


R.(N,T)  =  CB 


,a+l 


a  'X°ei-lTi 

MMTj*p( - ^— ) 


6X0  N  m-1  ,  -\o0j  iT 

.  t  v  a  N't  o+i  i-i  m 

+  JiSA' d»)T»  P 


o+l 


0+1 


i  =  1,2, ...  ,N 


(2.118) 


and  from  (2.111),  V.  is  calculated  to  be 


.o+l 


*o  Ti  a+i 

vi  *  - (ST 


(2.119) 


o+l 

Let  - — -  =  U°  A 

o+l 
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which  implies 


and  finally 


W.t. 

Vi  =  (-ei-lWi>~1  41  1exp(-U0,+1)Ua+1dU 


Wj  is  also  used  to  evaluate 
T, 

yi  *  4exP(--rtr — )dt 


which  can  be  shown  as 


W  T 

yi  =  ^o 1  W(-Ua+1)ciU 
From  (2,117)  and  (2.120)  we  have 

\o^ 

a+1  yi-l 


w: 


o+l 


o+l 


i  -  Wi-1 


A-o 

"o+r 


(8{ 


i-1 


-  0i-2>  = 


or 


r  0+1  _ 

i 


Wi-I1+ 


\0e 

a+T  yi-l 


W 


o+l 

1 


a+T 


i  =  2,3, ...,N 


(2.120) 


(2.121) 


(2.122) 


(2.123) 


(2.124) 
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The  following  computation  procedure  summarizes  the  results 

j|C 

obtained  in  this  section  to  evaluate  the  gradients  at  T  when  the  under¬ 
lying  failure  rate  has  the  form 

ri(t|  0i-i>  =  ei_1r(t) 

i-1 

and  the  age  factor  is  approximated  by  8.  ..  =  1  +  £  I  y. 

1-1  j=l  J 

STEP  0:  Set  i  =  1,  0O  =  1 

Ti- 

STEP  1:  Calculate  y.  =  SQl  F(t|8i_1)dt 

i 

STEP  2:  Compute  6.  =  1+  Z  y.  =  0.  .+  y. 

l  j=1  i  i-i  1 

STEP  4:  If  i  <  N,  set  i  =  i+1  and  go  to  Step  1;  else  go  to  Step  5 

STEP  5:  Calculate  matrix  d^;  i,j  =  1,2,  ...,N  using  eq.  (2.112) 

STEP  6:  Calculate  R^N.T);  i  =  1,2,...,N,  from  equations  (2.107) 

and  (2.11)  using  the  results  obtained  in  Step  5 

STEP  7:  Compute  L.(N,T);  i  =  1,2,...,N  from  equation  (2.114) 

STEP  8:  Calculate  R(N,T)  using  equations  (2.92),  (2.110)  and  the 

results  obtained  in  Step  6 

N 

STEP  9:  Calculate  L(N,T)  =  1  y. 

i=l  1 

STEP  10:  Evaluate  C.(N,T)  using  equation  (2.115)  and  the  results 
obtained  in  Steps  6  through  9 

Note  that  if  the  underlying  lifetime  distribution  is  Weibull,  then 
the  following  changes  are  made  to  reduce  computation  time. 

1.  W.  is  calculated  from  either  (2.120)  or  (2.124)  before  Step  2  is 
evaluated,  and  it  is  used  to  compute  V.  from  equation  (2.121)  and  y. 
from  equation  (2.123). 

2.  R^N,^)  is  evaluated  by  using  equation  (2.118). 
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Generally  any  optimum  search  technique,  when  N=l,  requires  an 
initial  interval  of  uncertainty  (X^.Xr).  We  may  define  the  lower  limit 
by  a  time  before  which  there  would  be  a  small  probability  that  the 
system  could  have  failed,  i.e. 


Pr{Tx  £  a}  £  6 

For  the  Weibull  distribution,  we  have 


A0a 


a+1 


a+ 1 


S  1-6 


or  finally 
a  £ 


-(g+l)ln(l-6) 

Ao 


1 

a+l 


=  XT 


With  a  similar  argument,  we  find  that 


XR  £ 


-(g+l)ln6 

Ao 


1 

a+1 


-  Xr 


(2.125) 


(2.126) 


(2.127) 


For  6  =  0.001,  equations  (2.126)  and  (2.127)  become 


0.001(q+l) 

Ao 


1 

a+1 


(2.128) 


1 

and  XR  = 

If  the  optimal  solution  converges  to  XR,  then  from  Figure  2.2(a) 
we  can  suspect  that  optimal  Tx(l)  is  infinity. 


6.9(q+l) 

Ao 


a+1 


(2.129) 


Numerical  Example  2.4 

We  now  consider  the  numerical  example  2.3,  but  we  employ  the 
approximation  method  to  solve  it.  The  results  are  shown  in  Figure  2.5. 
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_# 

The  optimal  policy  I  in  this  case  is:  N*  =  6,7  with  C(N,T  )  =  15.49  and 
T*  =  (.96  .90  .85  .81  .78  .74)  or  T*  =  (.94  .89  .83  .78  .75  .72  .69). 
We  have  used  the  conjugate  gradient  technique  and  the  Davidson, 
Fletcher,  and  Powell  method.  Obviously,  the  conjugate  gradient  tech¬ 
nique  needs  less  computation  time.  The  result  of  lemma  2.1  is  also 

applied  to  this  case,  and  computations  could  have  been  reduced  by 

*  * 

obtaining  an  estimate  for  N  (say  N'  ).  This  is  achieved  by  finding 
the  minimum  of  C(N,»). 


C(N.oq) _ 

N  Simulation  Approx. 


N 


1 

2 

3 

4 

5 

6 

7 

8 
9 


21.45 

18.63 

17.63 
17.54 
17.62 
17.99 
18.23 
18.49 
18.88 


21.45 

18.53 

17.91 

17.84 

17.96 
18.17 
18.43 
18.69 

18.97 


Table  2.2.  Comparison  between  simulation  and 
approximation  when  T  -*•  » 


Table  2.2  shows  that  N'  =4,  which  could  have  served  as  an  initial 
value  for  N  and  iterations  1,  2,  and  3,  could  have  been  avoided.  The 
expected  value  of  each  period  in  this  case  is 


a  1 


i-1 

where  0.  *  =  1+  6  i  p.. 

1  j=l  J 

Table  2.2  compares  the  limiting  behavior  of  C(N,®),  as  obtained  by 
simulation,  to  the  values  that  result  from  the  "mean  value”  approxima¬ 
tion  (numerical  examples  2.3  and  2.4). 


2.10  AVAILABILITY 


Throughout  our  discussion  about  policy  I,  we  assumed  that  the 
durations  of  times  for  maintenance  and  renewal  are  negligible.  In  a 
sense,  this  does  not  seem  to  be  justifiable.  But  by  interpreting  the 
costs  Cg,  CQ  and  Cg  as  the  expected  replacement,  maintenance,  and 
breakdown  times,  we  may  formulate  the  type  I  policy  in  terms  of  the 
steady  state  availability  (limiting  efficiency)  A(N,T). 


A(N,T)  = 


N  T, 

.f/o  FiCt|e.-1)dt 


N 


N  T. 


CRHN-l)Co+CBliFi(Tilei.i)+ii;o1Fi(t|ei.i)dt 


[1  +  C(N  ,T  )]-1 


Therefore,  maximizing  A(N,T)  is  equivalent  to  minimizing  C(N,T). 


CHAPTER  3 


POLICY  II 

SCHEDULING  FIXED  TIMES 

FOR  MAJOR  REPAIRS  WHEN 
INTERMEDIATE  MINOR  REPAIRS 

ARE  ALLOWED 


3.1  INTRODUCTION 

This  chapter  is  concerned  with  policy  II  and  its  optimalities  when 
the  major  repairs  are  scheduled  and  minor  repairs  are  carried  out  when 
failures  occur.  Minimizing  the  long-term  average  cost  per  unit  time  is 
used  as  a  criterion  to  optimize  this  strategy.  A  general  definition  of 
the  type  II  policy  is  found  in  Section  3.2.  Two  different  failure  rates 
(A  and  B  types)  are  defined  in  Section  3.3.  The  problem  is  formulated 
in  terms  of  nonlinear  and  dynamic  programming  in  Sections  3.5  and  3.6 
respectively.  Necessary  and  sufficient  conditions  of  optimality  for  this 
problem  are  derived  in  Section  3.7.  Section  3.8  contains  a  review  of 
the  state-dependent  model.  Sections  9-11  are  devoted  to  various  age- 
dependent  models  and  their  optimalities.  A  steady  state  availability 
formulation  is  found  in  Section  3.12. 

3.2  POLICY  II 

Unlike  policy  I,  which  is  suitable  for  simple  equipment  as  des¬ 
cribed  in  Chapter  2,  policy  II  is  applied  to  expensive,  complex,  and 
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multi-item  systems  (e.g.,  computers,  airplanes,  bulbs  in  a  factory, 
etc.)-  In  maintaining  a  large  and  complex  system  over  an  infinite  time 
span,  the  failure  of  a  single  component  unit  does  not  necessarily  call 
for  replacing  the  entire  system.  Instead,  the  system  can  be  restored 
to  operation  by  replacing  that  single  unit.  Generally,  any  maintenance 
action  of  this  type  is  called  "minimal  repair,"  after  which  the  system 
failure  rate  evolves  as  it  would  have  if  no  failure  had  taken  place.  A 
system  having  an  increasing  failure  rate,  after  a  certain  amount  of 
operation  time,  is  subject  to  renewal  and/or  overhaul  (major  repair). 
Our  objective  is  to  determine  the  optimum  maintenance  time  intervals 
(time  between  two  successive  major  repairs)  and  renewal  cycle  or  simply 
a  cycle  which  is  defined  to  be  the  time  between  two  successive  replace¬ 
ments.  For  a  complex  system  and  when  the  first  major  repair  is  a 
replacement,  a  strategy  known  as  policy  II  has  received  the  most 
attention  in  the  literature  [Barlow  and  Proschan  (1965)]. 

The  following  is  a  generalization  of  the  type  II  policy: 

Policy  II.  The  system  is  repaired  for  the  i1^1  time  at  age 
i 

t.  =  l  T.,  provided  i  <  N;  it  is  replaced  when  i  =  N.  In  the 
1  j=l  J 

case  of  failure  between  maintenance  actions,  a  minimal  repair 

is  carried  out.  This  process  will  be  repeated  indefinitely. 

Our  assumptions  made  in  Section  2.1  regarding  zero  time  duration 
for  replacements,  major  and  minor  (minimal)  repairs,  lack  of  queuing 
problem  for  maintenance,  system's  state  being  observable,  and  having 
increasing  failure  rates  are  also  assumed  for  policy  II. 

This  strategy  considers  two  kinds  of  repairs: 


1.  In  the  case  of  failure,  the  system  can  be  restored  to  operation 
by  minimal  repairs,  after  which  the  failure  rate  evolves  as  it  would 
have  if  no  failure  had  taken  place. 

2.  Planned  major  repairs  in  TltT2, . . .  ,T^_^  units  of  time  from 
the  most  recent  corrective  operations  (major  repairs)  and  then  replace¬ 
ment  at  age  t^,  where  t.  is  defined  to  be 

tj  =  2  T.  (3.1) 

1  j=l  J 

The  maintenance  costs  in  this  model  are  defined  as  CR  for  each 
replacement,  which  is  performed  at  age  t^j,  Cj^  for  each  major  repair  at 
ti/t2,  •  •  •  ,tN-1  and  CQ  for  each  minimal  repair  in  time  interval 
i  =  1,2,...,N.  These  costs  are  assumed  to  be  constant  despite  the  fact 
that  they  might  well  be  functions  of  age,  number  of  repairs  (major 
and/or  minor),  and  many  other  factors. 

3.3  MODELS  A  AND  B 

After  each  major  repair  that  is  carried  out  to  reduce  the  hazard 
rate  (and  hence  the  number  of  minimal  repairs),  the  system  does  not 
completely  recover,  and  breakdowns  are  more  likely  to  occur  than  if  a 
replacement  had  been  performed.  Mathematically,  this  can  be  described 
by  a  failure  rate  function  after  the  (i-l)st  major  repair,  r^t), 
satisfying 

ri(ti_l+At)  S  ri.1(ti.2+At)  2  ...  S  r1(t0+At)  (3.2) 

The  function  r.  might  depend  explicitly  on  i-1,  V^t^ , ...,t1. 
Figure  3.1  is  a  graphical  representation  of  relation  (3.2).  In  this  case, 
major  repairs  will  reduce  the  failure  rates  such  that 
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ri(Vi)  =  ri_i(ti_2)  =  •"  =  ri^o)  (3.3) 

We  refer  to  this  model  as  "Model  A."  Another  illustration  is  shown 
in  Figure  2.1.  When  relation  (3.3)  is  ignored,  one  can  define  "Model  B" 
for  which  a  typical  example  is  represented  in  Figure  (3.2). 

3.4  PROBLEM  FORMULATION 

Let  0.  =  gi(T1/T2,...,T._1)  be  a  positive  real  function  to  describe 

of 

the  past  history  of  the  system  immediately  after  the  (i-1)  preventive 
maintenance  (major  repair)  action.  This  definition  will  permit  us  to 
formulate  the  type  II  policy  when  the  failure  rate  depends  on  the  age  of 
the  system  and/or  number  of  repairs. 

The  failure  characteristic  of  the  system  can  be  expressed  by  the 
conditional  failure  distribution  function  Fi(t|0i  l),  which  is  related  to 
the  conditional  hazard  rate  r.(t|6i  l)  and  FA(t | ©j  .^)  is  assumed  to  be 
strictly  increasing  in  Tx,T2, . . . /Tj^,  t,  and  possibly  (i-1).  We  assume 
that  Fj(t|0._1)  and  r^tie^)  are  both  continuous  and  differentiable. 

The  expected  total  cost  per  renewal  cycle,  denoted  by  R(N,T),  is 
computed  by  adding  the  replacement  cost,  (N-l)  major  repair  costs,  and 
the  expected  minimal  repair  cost  caused  by  failures.  Based  on  our 
notations  defined  earlier,  we  have 

N 

R(N,T)  =  CR  +  (N-1)Cq  +  CM  Z^E[N(t._1,ti)]  (3.4) 

where  N(t._1,tj)  is  defined  to  be  the  number  of  failures  (and  hence  the 
number  of  minimal  repairs)  occurring  in  (t._^,tp,  and 

T  *  . . v 


(3.5) 
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By  definition,  the  probability  of  a  failure  occurring  in  (t,t+dt)  is 
rj(t|0._i)  +  0(t).  Barlow  and  Hunter  (1961)  and  Barlow  and  Proschan 
(1965)  have  shown  that  the  expected  number  of  failures  in  (^.pt.)  can 
be  expressed  as 

ElNCtj.j.t,)]  =  r|(t-tj_1|01-1)dt 

=  Jo1  ti(x|e._1)dx 

-  ziCTilei_1)  (3.6) 

This,  of  course,  follows  from  the  fact  that  the  minimal  repairs 
during  time  interval  (t._^,t.)  do  not  disturb  the  system's  failure  rate. 
In  other  words,  after  a  minimal  repair,  we  assume  that  the  failure  has 
never  taken  place. 

Thus,  by  using  equation  (3.6),  we  can  rewrite  the  identity  (3.4) 
as 


N  T. 


R(N,T )  =  CR  +  (N-1)C0  +  C^j  I  S0  ri(t|6..1)dt 


(3.7) 


The  expected  length  of  a  cycle  is  the  total  age  of  the  unit 

N 

L(N,T)  =  tN  =  IT,  (3.8) 

w  i=l  1 

The  expected  long-term  average  cost  for  this  model  is  just  the 
expected  cost  during  a  cycle  divided  by  the  expected  length  of  a  cycle. 
Therefore,  the  expected  cost  per  unit  time  is  obtained  by 


C(N,T) 


_R(N,T) 

L(N,T) 


CrHN-DVCm  I=1  s} 

- FT - 

I  T. 


ri(t,ei-l)dt 


(3.9) 
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Note  that  the  hazard  rate  is  affected  by  age  and/or  repair  only 
after  each  major  repair,  and  it  is  not  disturbed  by  minor  repairs  which 
take  place  before  the  next  major  repair. 

We  now  turn  our  attention  to  the  optimization  of  preventive  main¬ 
tenance  policy  II. 

3.5  OPTIMAL  POLICY  II 

In  order  to  find  an  optimal  policy  II,  we  seek  a  set  of  positive 
time  intervals  (T.;  i  =  1,2,...,N),  which  minimizes  the  expected  long¬ 
term  average  cost  per  unit  time  given  by  (3.9).  Mathematically,  this  is 
expressed  by  the  following  nonlinear  programming  problem  (NPII). 

Find  N*  and  T*  =  (T1(N),T2(N) . TN(N» 


To  minimize: 

C(N,T)  =  (3.9) 

L(N,T) 

Subject  to: 

T.  £  0  i  =  1,2,. ..,N  (3.10) 

NU  (3.11) 

N  Integer  (3.12) 


In  a  way  similar  to  (NPI)  given  in  Chapter  2,  we  have  three 
special  cases  to  consider. 

Case  I.  No  major  repairs  (N=l) 

When  major  repairs  do  not  result  in  any  improvement  or  the  only 
possible  repair  is  the  minimal  repair,  then  N*  =  1.  Barlow  and  Hunter 
(1960)  and  Barlow  and  Proschan  have  studied  this  model.  They  intro- 
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duced  this  notion  as  "periodic  replacement  with  minimal  repair  at  failure." 
For  N  =  1  equation  (3.9)  yields 


C(1,T) 


°R  +  °M  ri^dt 
Ti 


(3.13) 


Setting  its  first  derivative  to  2ero,  we  get  the  first  necessary 
condition  for  a  minimum 


CMTlr1(T1)  -  [CR*CUSQ  rx(t)dt]  =  0  (3.14) 

or 

A  Ti  CR 

Gi(Tj)  t  Tjr^Tt)  -  J0  rx(t)dt  =  ^  (3.15) 

Therefore,  the  optimum  time  interval  between  replacements  is 
obtained  by  Tx(l),  which  is  the  solution  to  equation  (3.15).  The 
minimum  cost  is 

C*(1,T1(1))  =  CMr1(T1(l))  (3.16) 

Barlow  and  Proschan  (1965)  have  sh  wn  that  Tx(l)  exists  and  it  is 
unique.  This  can  be  seen  by  differentiating  the  left-hand  side  of 
(3.15)  to  get 


dG^TO 

“377“  =  TWi* 


(3.17) 


which  is  positive  when  r(t)  is  increasing  in  t.  The  right-hand  side  of 
(3.15)  is  constant  and  hence  Tx(l)  is  unique.  Figure  3.3  illustrates  a 
typical  expected  cost  function,  CQ.TO. 

As  an  example,  suppose  the  failure  rate  has  the  Weibull  form 


ri(t)  =  \Qta 


t  >  0  ;  a  >  0 


(3.18) 


Equation  (3.15)  can  be  written  as 


a+1 

1 


\>  rp  a+1  _ 

^rTi 


°M 


or 


V 

a+1 


»,a+l 
1  i 


(3.19) 


(3.20) 


and  the  optimum  time  period  between  replacements  is  given  by 


Ti(l)  = 


(a+l)CR 


I  a+1 


aAoSl 


(3.21) 


The  resulting  minimum  expected  cost  is  obtained  by  using  identity  (2.16) 


C*(l,  Tid))  =  C^T^l))"  = 


M  O 


(a+l)CR 


1 

la+1 


(3.22) 


Case  II.  Rapid  Aging 

This  case  is  similar  to  Case  II  described  in  Chapter  2.  It  is 

applied  when  major  repairs  have  a  destructive  nature  or  the  system 

deteriorates  very  rapidly.  As  mentioned  in  Chapter  2,  the  optimal 

* 

number  of  major  repairs  is  N  =1.  Note  that,  in  this  case,  minimal 
repairs  are  performed  for  continuing  operations,  and  they  do  not 
disturb  the  failure  rate. 

Case  III.  Slow  Aging 

Among  the  two  kinds  of  repairs,  minor  repairs  do  not  affect  the 
system's  failure  rate,  but  major  repairs  will  result  in  improvement.  This 
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case  assumes  that  after  each  major  repair,  the  system  becomes  as  good 
as  new.  A  mathematical  representation  of  this  model  is  described  in 
terms  of  hazard  rate 

ri(T|0i-l)  =  r(t)  (3-23) 

and  the  resulting  long-term  average  cost  per  unit  time  given  by  relation 
(3.9)  becomes 


C(N,T) 


cr+(n‘1>co+cm  !/Iir(t>dt 


N 

I 

i=l 


(3.24) 


Setting  its  first  derivatives  to  2ero  to  get  the  first  condition  of 
optimality 


N 


N  T. 


CMr<Tj>  ,3/1  -  (CrHN-1)C0*Cm  Z^o'rCDdt]  =  0 
j  —  1,2,.. .,N 


N  N  T. 


or  r(T.)  I  T.  -  Z  f  lr(t)dt  = 
J  i=l  1  i=l  0 


.  V^^o 


°M 


j  “  1,2,. ..,N 


From  (3.25),  minimum  cost  is  found  to  be 

C(N,T*)  =  CMr(Tj(N))  j  =  1,2 . N 

which  yields 

r<Tj(N))  =  r(Tj(N))  i,j  =  1,2 . N 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


Since  r(t)  is  an  increasing  function  of  t,  dearly 


T.(N)  =  T.(N)  =  T  i,j  =  1,2 . N 

Combining  (3.29)  und  (3.26),  we  get 


(3.29) 


<t>  CR+(N-1)C 

NTr(t)  -  J*  r(t)dt  =  - 5 


(3.30) 


Thus  the  optimum  time  interval  between  two  major  repairs  (or  possibly 
replacements)  is  obtained  by  T,  which  is  the  solution  to 


T  Cr+(N-1)C 

G(T)  £  Tr(T)-  r(t)dt  =  - 2  = 

M 


Differentiating  G(T)  with  respect  to  T  yields 


h(N) 


(3.31) 


dG 


=  Tr'(T)  2  0 


(3.32) 


The  right-hand  side  of  (3.31),  h(N),  is  increasing  in  N  if 

CR  >  CQ  (3.33) 

In  an  argument  similar  to  the  one  given  for  policy  I  (Chapter  2,  special 
Case  III),  we  conclude  the  following: 

1.  If  major  repair  cost  C  is  more  than  replacement  CD,  then 

O  I\ 

* 

the  optimum  number  of  preventive  maintenance  actions  is  N  =1  (the 
same  as  Case  I). 

2.  If  major  repair  cost  is  less  than  replacement  cost,  then  N*  =  », 
and  the  optimum  time  interval  between  major  repairs  is  T,  which  is  the 
solution  to 

C 

Tr(T)  -  r(t)dt  =  ^  (3.34) 
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3.  The  optimum  number  of  planned  major  repairs  is  independent 
of  minimal  repair  costs,  C^.  But  larger  time  interval  for  each  scheduled 
major  repair  is  obtained  by  smaller  C^. 

4.  When  CR  =  CQ,  then  N  can  be  any  positive  integer  number. 

The  above  results  are  somehow  obvious.  But  for  state,  age, 
and  state-age-dependent  policy  II,  we  assume  that  Cr>  >  C  .  This 
allows  us  to  investigate  the  existence  of  nontrivial  optimal  solutions. 
Case  III  is  a  special  case  of  state-dependent  model,  developed  by 
Nguyen  and  Murthy  (1981),  which  we  shall  briefly  discuss  in  Section 
3.7. 

Before  considering  some  more  general  models,  a  reformulation  of 
the  optimal  policy  II,  using  the  dynamic  programming  technique,  is 
presented. 

3.6  DYNAMIC  PROGRAMMING  FORMULATION 

This  problem  is  analogous  to  the  one  discussed  in  Section  2.4. 
Equations  (2.39),  (2.40),  (2.41),  (2.45),  and  (2.46)  still  hold.  Identity 
(2.42)  will  change  to 

T 

CK  =  a-r'K>CR1-nKVCM  4  r<tleK-l>dt  -  PTK  «.35) 

where 

nK=l  if  1  £  K  <  N 

=  0  if  K  =  N 

A  typical  transition  equation  can  be  defined  by  0^  when  it  is  pro¬ 
portional  to  the  system's  age  at  the  time  of  the  most  recent  major  repair 

0K  =  8K-1  +  TK  k  *  0 

eo  =  1  (3.36) 
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where  €  is  a  constant  to  determine  how  fast  the  system  deteriorates. 
This  relation  is  especially  suitable  when  hazard  rate  has  the  form 

ri(tlei-l)  =  e^rCt)  (3.37) 

i  =  1,2, . . . ,N 

which  reflects  the  A-type  model  defined  in  Section  3.2.  For  B-type 
model  eQ  can  take  on  any  nonnegative  value. 

When  failure  times  follow  the  Weibull  distribution,  obtaining  optimal 
policy  II  by  means  of  dynamic  programming  requires  numerical  compu¬ 
tations.  At  each  stage,  an  equation  involving  parameters  like  6^  and  0 
must  be  solved.  Bounds  on  0^  and  0  may  reduce  inefficiency.  A  more 
efficient  method  is  to  apply  a  search  technique  (e.g.  gradient  method). 
For  this  to  be  the  case,  the  nature  of  optimal  solution(s)  must  be 
known. 


3.7  EXISTENCE  OF  MINIMUM 

We  have  seen  that  there  always  exists  a  minimum  when  N  =  1. 
The  following  notations  are  used  to  investigate  the  existence  and  the 
nature  of  any  stationary  point  for  (NPII)  when  N  £  2. 

Let  the  first  and  second  partial  derivatives  of  total  cost  per  cycle, 
R(N,T),  be 


r.(n,¥)  = 

and  RjjCN.T)  = 

Also  let 

=  ac(N,T) 


i  =  1,2, ... ,N  (3.38) 

i,j  =  1.2 . N  (3.39) 


C^N.T) 


i  =  1,2 . N 


(3.40) 
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and  C.j(N,T) 


92C(N,T) 

aT^. 


i,j  =  1.2 . N  (3.41) 


be  the  first  and  second  partial  derivatives  of  C(N,T). 

For  rK^leK-l^'  s^nce  t^le  subscript  represents  the  state  of  the 
system  (number  of  major  repairs  performed),  we  employ  superscripts  to 
distinguish  the  first  and  second  partial  derivatives.  Thus 


m  8rK<t,eK-l> 

■kWr-i)  *  KarK1 


i  =  1,2,...,K-1 


(3.42) 


rK  }  «I9K- 


32rK:(t|0K:_1 ) 

i)  =  — tr.aT'""  =  1>2 . K_] 


(3.43) 


Note  that  rK  is  independent  of  TV,Tj£  , , . . .  ,T^. 

To  minimize  C(N,T),  we  set  its  derivatives  equal  to  zero 


Cj(N,T)  = 


R.(N,T)  R(N,T) 


N 


j=i 


N 

(  l  T.)2 
j=l  J 


=  0 


(3.44) 


i  =  1,2,. ..,N 


The  first  necessary  condition  for  minimum  is  obtained  by 
N 

R.(N,T)  XT.  -  R(N,T)  =  0  (3.45) 

1  1=1  ] 

or 

R.(N,f)  =  =  C(N,T)  (3.46) 

X  T. 

j=l  ] 

i  =  1,2 . N 

which  gives 

R.(N,T)  =  R.(N,T)  i,j  =  1,2 . N  (3.47) 
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Therefore,  the  optimal  solution  is  such  that  the  change  in  total  expected 
minimal  cost  per  cycle  will  remain  the  same  when  both  T.  and  T.  (tfj) 
are  incremented  by  one  unit.  From  (3.7)  R.(N,T)  is  found  to  be 


Rj(N,T)  = 


Note  that 


N 

2 

j=i+l 


( 


M 


N  T. 


ri<Ti'6M>  +  j4+/«S(,>(t|eH)dt 


i  =  1,2,. ..,N 


)  =  0  when  i  =  N. 


(3.48) 


The  minimum  long-term  average  cost  per  unit  time  can  be 
expressed  in  terms  of  the  failure  rate  at  the  end  of  the  last  stage  by 
using  relations  (3.48)  and  (3.46)  when  i  =  N. 

C(N,T*)  =  CMrN(T^|eN-1)  (3.49) 


One  can  substitute  (3.48)  and  (3.7)  into  equation  (3.45)  to  get 
N 


ri(Ti'9i-i)  *  rjl)  «iVi)dt 


N 

2  T 
j=l  1 


-  CR  -  (N-1)C0  -  Cy  rj(t|9j_1)dt  =  0 


(3.50) 


which  can  be  simplified  to 


N  T. 


N 

2  T. 
j=l  3 


N  T. 


-  .J/o  rjCtiej-iXtt  = 


_  Cr+(N-1)Cc 


°M 


i  =  1,2 . N 


(3.51) 
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Setting  i  =  N  yields 
N 


N  T. 


. ,  Cr+(N-1)C 

rN<TNl0N-l>^Ti  -  jf/o'  'j(tlVl)dt  — 


(3.52) 


Combining  (3.47)  and  (3.48)  leads  to 


ri<Ti'0i-l>  -  rj<Tjl6j-l>  = 

2  4KrK)(t|6K-l)dt  ‘  *  4K  rK)(t|6K-l)dt 

K=j+1  0  K  K  1  K=i+1 0  K  K  1 


(3.53) 


i,j  =  1,2,.  ..,N 


(3.53)  along  with  (3.52)  constitutes  a  system  of  N  simultaneous,  non¬ 
linear  algebraic  equations  with  N  unknowns  which  can  be  solved  to 
obtain  the  optimal  solution(s).  This  is  also  achieved  by  solving  (3.51) 
but  with  more  calculations. 

To  evaluate  the  second  partial  derivatives,  we  define  K  =  Min(i,j) 
and  m  =  Max(i,j)  for  each  pair  of  (i,j)  to  obtain 


R^N.T) 


N 

again  2  (  ) 
n=m+l 


=  C 


M 


NT.. 

^riK)(Tmlem  i)  +  2  Snn  r(lj)(t|e„  ,)dt 

m  m  m-1'  n=in+l  0  n-1 


k,j  =  1,2 . N 


0  when  m  =  N.  C..(N,T)  is  found  by  using  eq. 


(3.54) 

(3.44). 


N 


Cij(N'T)  = 


R..(N,T)  2  T  +  R.(N,T)  -  R.(N,T) 
m=l m  1  ] 


N 

(  2  T  )2 
m=l  m 


N 


2[R.(N,T)  2  Tm  -  R(N,T)] 
m=l 


N 

(  2  T  )» 
m=l  m 


(3.55) 


i,j  =  1,2,... ,N 
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Sufficient  conditions  for  minimum  are  considerably  complex  if  we 
deal  with  Hessian  matrix 


H  =  [C..(N,T>] 

where  C.j(N,T)  is  given  by  (3.55).  But  these  conditions  can  be  sim¬ 
plified  if  we  evaluate  H  at  a  stationary  point.  Since  (3.45)  and  (3.47) 
will  hold  for  the  optimal  solution,  we  can  substitute  them  into  the 
equation  (3.55)  to  get 


C..(N,T*) 


R..(N,T  ) 
N 

I  T  (N) 
m=l  m 


(3.56) 


The  following  analysis  is  used  to  discuss  the  sufficient  conditions, 
for  which  the  N-dimensional  function  C(N,T)  is  optimal  [Kiepert  (1910)]. 
Let  the  set  of  determinants  |D.|;  i  =  1,2,..., N  be 


ID, 


* 

Rxi 


* 

R2i 


* 

Ril 


* 

R12 


* 

R22 


*12 

L 


R 


li 

L 

* 

L2i 

"1 

L 


* 

Rii 


(3.57) 


Then  we  have  the  following 


1.  |Dj|  <  0  for  i  =  1,3,...  and  |D.|  >0  for  i  =  2,4,...  indicates 
the  presence  of  a  relative  maximum  at  T*. 

$ 

2.  |Dj}  >  0  for  i  =  1,2,...,N  represents  a  minimum  at  T  . 
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3.  Failure  to  satisfy  either  of  these  two  conditions  indicates  the 
presence  of  saddle  point. 

Due  to  the  differentiability  of  R(N,T)  and  C(N,T)  with  respect  to 
TX,T2 . Tn,  we  have 


Cj.(N,T)  =  Cji(N  ,T ) 

(3.58) 

and  RjjCN/T)  =  R^N,?) 

(3.59) 

Define 

|D!|  =  aVlD.I 

i  =  1,2 . N 

(3.60) 

* 

Since  L  >  0,  we  can  simplify  the  second  order  optimality  condition 
by  means  of  |D||.  In  other  words,  any  relation  that  satisfies  |D.|  will 
also  satisfy  |  Dj  | . 


Rxi  Rxa 

|Dj|  =  R2i  R22 

•  • 

•  • 

•  • 

*  * 

Ril  Ri2 

Theorem  3.1. 

Under  the  folkring  conditions: 

HI:  r.(t|6._p  is  strictly  increasing  in  T1#T2,. . .  and  t 
H2:  ^ C 1 1 ©m- 1  ^  ®  *or  a*  i  and  m 

C(N,f )  has  a  positive  relative  minimum  with  respect  to  Tj  in  the 
interval  (0,»),  and  the  solution  to  the  first  necessary  condition 
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R.(N,T)  =  C(N,T)  (3.46) 

cannot  be  a  relative  maximum. 

Proof.  H2,  (3.54),  and  (3.61)  imply  that  |Dil  >  0,  and  hence  the 
relative  maximum  point  cannot  exist.  C(N,T)  is  a  continuous  function 
of  Tt;  i  =  1,2,...,N;  and  it  approaches  infinity  as  T  ■*  0  or  ■*  ».  Note 
that  C(N,T)  may  possess  a  saddle  point. 

The  preceding  theorem  allows  us  to  employ  various  optimum  search 
techniques  to  find  the  minimum.  The  next  theorem  is  useful  for  finding 
proper  initial  values  for  T1#T2, . . .  ,T^. 


Theorem  3.2 

If  (i)  r.(t|0._1)  is  increasing  in  TlfT2,...,TN  and  t,  and 
(ii)  r(im)(t|0i_1)  =  r<n>(t|0H)  for  m  =  l,2,...,i-l  and  n  =  1,2,..., J-l, 
then  failure  rates  evaluated  at  optimal  T.  are  ordered  according  to 


ri<Ti'W  * 


i  >  J 


(3.62) 


Proof.  Without  loss  of  generality,  we  let  j  =  i-1,  and  from  (3.53) 
and  hypothesis  (ii),  we  get 


(3.63) 


Condition  (i)  implies  that  the  right-hand  side  of  (3.63)  is  positive 
and  hence 


ri(Ti|6i-l)  =  ri-l(Ti-ll0i-2)  1  =  l'2 . N 

As  an  example  for  which  conditions  of  Theorem  3.2  satisfy,  con¬ 
sider  the  Weibull  failure  rate 
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ri<t»ei-i>  =  Vi-/  “  * 1 

i-1 

®i_l  =  1  +  €  I  T. 

11  j=l  ] 

Note  that  =  /.n)  =  6t“;  a  2  1. 

So  far,  we  have  studied  the  general  policy  II.  Now  we  turn  out 
attention  to  state-  and  age-dependent  type  II  policy  separately. 


3.8  STATE-DEPENDENT  MODEL 

In  this  case  of  policy  II,  the  state  represents  the  number  of  major 
repairs  performed.  The  failure  rate  at  each  stage  is  an  increasing 
function  of  the  number  of  major  repairs  already  carried  out  and  the 
time  from  the  most  recent  major  repair.  Mathematically 

ri(t|6i-l)  =  ri(t)  t  >  0;  i  =  1,2 . N  (3.64) 

From  (3.9),  the  long-term  average  cost  per  unit  time  is  found  to  be 

Or^N-DGo-Om  X/o  V>dt 

C(N,T )  =  - -  (3.65) 

■  XTi 
1=1  1 

D.G.  Nguyen  and  D.N.P.  Murthy  (1981)  have  studied  this  model  for 
both  policy  I  and  policy  II.  They  have  shown  that  these  two  models 
have  similar  behaviors.  Their  results  are  summarized  as  follows: 


1.  Differentiating  (3.65)  with  respect  to  t.  and  equating  it  to  zero 
shows  that  the  optimum  preventive  schedule  times  satisfy 


r.(T.)  =  ri(Tx) 


i  =  2,3,.. .,N 


N  NT, 

and  ri(Tt)  IT,  -  I  J  Wodt 
i=l  1  i=l  0  1 


Cr+(n-dc0 

°M 


(3.66) 


(3.67) 
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and  the  resulting  minimum  cost  is 

C(N,T*)  =  C^r-CT.)  (3.68) 

2.  Assumptions  (i),  (ii),  and  (iii)  made  in  Section  2.5  will  guar¬ 
antee  the  existence  and  uniqueness  of  an  optimal  N  and  optimum  planned 
major  repair  time  intervals  T.(N);  i  =  1,2,..., N.  For  a  fixed  N,  T\(N) 
has  been  shown  to  be  decreasing  in  i. 

3.  When  failure  times  are  assumed  to  be  given  by  Weibull  distri¬ 
butions 

fj(t)  =  \jt“exp(-  ta+1)  a  i  1  ;  t  >  0  (3.69) 


Then  T^(N)  can  be  found  analytically  by  letting 


Ax  \l/cf 


ni  =  l  k: 


which  yields 
TX(N)  = 

and 


(n-dCq+Cr 
N 

L  c^xMi^D. 


-il/(a+l) 


T.(N)  - 


niTi(N); 


i  =  1,2 . N 


(3.70) 


(3.71) 


4.  Nguyen  and  Murthy  have  also  shown  that  the  computation 
algorithm  used  for  policy  I  (see  Section  2.5)  can  be  applied  to  policy  II, 
except  that  here  Step  (ii)  is  computed  by  solving  (3.66)  and  (3.67), 
and  Step  (v)  is  obtained  by  (3.65).  But  our  heuristic  procedure  to  find 
approximate  N  for  policy  I  cannot  be  used  for  policy  II. 
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3.9  AGE -DEPENDENT  MODEL 

Unlike  policy  I  in  which  the  age  at  the  time  of  repair  was  a 
random  variable,  here  the  age  of  the  system  is  known  and  constant 
when  a  major  repair  is  performed,  and  it  can  be  represented  by 

i 

t.  =  IT.  i  =  1,2 . N  (3.72) 

1  j=l  ] 

where  tN  indicates  the  replacement  time.  We  consider  two  types  of 
age-dependent  policy  II  as  follows: 

1.  One  corresponds  to  "Model  A"  defined  in  Section  3.2.  In  this 
case,  the  hazard  rate  at  the  beginning  of  each  stage  is  given  by 

ri(t|0i-l>  =  ei-lr(t)  1  =  1'2 . N  (3‘73) 

t  £  0;  60  £  1 

As  mentioned  earlier,  0._^  is  an  increasing  function  of  T1/T2, •  •  • 

2.  In  the  second  model,  the  failure  rates  have  the  form 

ri(t|0i-i>  =  r(t)  +  0^  i  =  1,2 . N  (3.74) 

t  §0  ;  0o^O 

which  is  a  particular  case  of  "Model  B." 

3.10  MODEL  A 

The  age  factor  for  this  case  may  be  defined  as 

i 

0.  =  1  +  t.  =  1  +  €  Z  T.  (3.75) 

1  1  j=l  ] 

where  €  is  a  positive  parameter  which  measures  the  degree  of  deteri¬ 
oration  as  the  system  gets  older.  Note  that  for  a  new  system,  6Q  =  1 
and  the  failure  rate  (3.73)  is  described  by  r(t).  From  (3.9),  we 
conclude  the  long-term  average  cost  for  unit  time 
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CR+(N-1>Co+CM 

lv\  -  _  1=1 


C(N,T)  = 


N 

•  Z  Ti 
i=l  1 


(3.76) 


The  derivative  of  r.  with  respect  to  Tj  in  this  case  is 


rii)(tl0j-1)  =  6r(t) 
=  0 


j  —  1/2,...  ,i-l 
j  =  i,i+l . N 


(3.77) 


which  can  be  substituted  into  relation  (3.51)  to  get  the  first  necessary 
condition  for  minimum 

NT.  N  N  T. 

[ri(Ti)6i-l  +  €j4+i/oJr(t)dt].ITj  -  __I  0j.i/o]r(t)dt 


i  =  1,2, ... ,N 


Equation  (3.53)  becomes 


T 

i+1. 


r(Ti+l)6i  *  r(Ti>ei-l  =  e/o  r(t)dt 


i  =  1,2, .. .  ,N-1 


(3.78) 


(3.79) 


As  an  example,  consider  a  system  whose  failures  follow  the  condi¬ 
tional  Weibull  distributions 


fi(t|0i-l)  =  Vi-/«P(-  ^  t°+1) 


(3.80) 


i  =  1,2,. ..,N 


0O  8  1;  t  i  0;  o  i  1 
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with  the  corresponding  conditional  hazard  rate 

ri<t|6i-l)  =  Vi-l1*  i  =  1.2 . N  (3.81) 

where  0.  is  given  by  (3.75).  In  this  case,  C(N,T)  has  the  following 
form: 


C(NfT)  = 


CR  +  (N-1)C0  + 


Vo 


N 


m-1 


<*+1  m-1 

m-1 


I  (1  +6  1 


N 

•  Z  Ti 
j=l  J 


Jzk 


T3>Tm+1 


(3.82) 


and  equation  (3.79)  will  change  to 


(1  +  €  I  T.  - 
j=l  ] 


a+l 


Ti+1>T?+1 


(1  + 


€  I  T.)Tf 
j=lJ  1 


(3.83) 


i  =  1,2, .. . ,N-1 

These  are  N-l  nonlinear  equations.  The  N^1  equation  is  the  result 
of  (3.78)  when  i=N: 


N-l 

(1+6  I  T.)T 
j=l  J 


nA 


Ti 


1 

a+l 


N  m-1 
I  (1+6Z  T.)T 
m=l  j=l  J 


a+l  _ 
m 


Cr+(N-1)Cc 

^oCM 


(3.84) 


Thus,  one  way  to  obtain  the  optimum  solution  is  to  solve  (3.83) 
and  (3.84)  simultaneously  for  TX,T2, . . .  ,TN.  But  according  to 
Theorem  3.1  saddle  points,  but  not  a  relative  maximum,  are  also  pos¬ 
sible  solutions .  The  following  theorem  may  help  us  to  distinguish 
minimum  point  from  saddle  points. 


Theorem  3.3 

If  the  underlying  distribution  is  Weibull  with  an  increasing  failure 
rate  and  scale  parameter 

\  *  V>t-i  *  V1*  «£Tj> 


(3.85) 
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Then  the  optimal  scheduled  repair  times,  T.,  for  policy  II  are  ordered 
as 

Ti  S  T2  ^  ...S  Tn  (3.86) 

Proof.  Consider  two  arbitrary  positive  numbers  a  and  b.  Let 


Ca  =  C(N,T)|T.  =  a,  T.+1  =  b 

and 

Cb  =  C(N,T)|T.  =  b,  T.+1  =  a 


where  each  Tj  in  Ca  and  Cb  are  respectively  the  same  for  j/i,  i+1.  It 
is  desired  to  show  that  for  any  a  >  b,  we  have 


Ca<Cb 


(3.87) 


N 


In  relation  (3.76),  we  can  see  that  the  total  age  l  T. 

j=l  J 

and  j>i+l)  remain  unchanged  after  interchanging  the  values 
T.+f.  Subtracting  from  yields 


Ca-Cb 


°MD 

N 

2  T. 
j=l  J 


and  0j(j<i 
of  T.  and 


where  D  is  defined  to  be 


D  =  r(t)dt  +  (9i_1+ea)J'o  r(t)dt 

‘  9i-i/or(t)dt  "  (ei.1+€b)4r(t)dt 

=  €la/Jr(t)dt  -  bjjr(t)dt] 


A 
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It  is  sufficient  to  prove  that  D  <  0.  Setting  r(t)  =  k  ta.  we  get 

L1  -  au  /V« 

1  "  5TT  (b  "a  > 


D=  eisrb' 


a  ua+l  b  q+1-j  _  ab  /Ka 


a+1 


but  we  assumed  that  a  >  b,  a  £  1  and  hence  D  <  0,  which  implies 
(3.87). 


Another  way  to  find  the  minimum  of  C(N,T)  is  to  apply  an 
optimum  search  technique.  This  in  turn  requires  analytical  expressions 
for  gradients  C.(N,T).  From  (3.48),  we  have 

*  NT. 

Ri(N.T)  =  CM[e.  r(t)  +  €  I  J  r(t)dt]  (3.88) 

j=i+l 

i  =  1,2 . N 


and  for  Weibull  case 


N 


¥N-T>  =  X0CM[(1*  e  i  T  )T“  *  Ij.  \  i^1] 

J=1  1  J=l+1  1 

i  =  1,2, ... ,N 

C.(N,T)  can  be  calculated  by 


(3.89) 


Cj(N.T)  = 


Rj(N,T)L(N,T)  -  R(N,T)L.(N,T) 


[l(n,T)]2 

i  *  1,2 . N 


(3.90) 


To  find  the  second  partial  derivatives  of  R(N,T),  let  m  =  max(i,j)  and 
hence 


RyCN.T)  =  eCMr<Tm)  ;  i/j 


*  cM9i-ir'<TP  :  *1 


(3.91) 


Note  that  conditions  HI  and  H2  given  in  Theorem  3.1  are  both 
satisfied. 

Numerical  Example  3.1 

Assume  that  the  failure  time  of  both  the  new  and  repaired  system 

have  Weibull  distributions  given  by  (3.80)  with  parameters  a=\Q=  =1. 

The  maintenance  costs  are:  0^=15. 0,  Co=1.0,  and  CM=0.3.  We  applied 

the  gradient  method  to  C(N,T)  (equation  (3.82)).  The  optimum  number 

* 

of  major  repairs  per  cycle  is  N  =8  with  the  minimum  expected  cost  per 
♦  * 

unit  time  of  C(N  ,T  )  =  2.88.  The  optimum  major  repair  times  are 

f*  =  (7.92,  0.88,  0.83,  0.80,  0.77,  0.74,  0.72,  0.70) 

_* 

Figure  3.4  illustrates  the  behavior  of  the  C(N,T  )  as  N  is  incre¬ 
mented.  It  is  a  convex  function  and  its  minimum  is  unique. 

* 

In  our  numerical  examples,  the  convexity  of  C(N,T  )  with  respect 
to  N  seems  to  be  a  sure  possibility.  In  addition  to  numerical  example  3.1, 
Figure  3.5  also  illustrates  this  property,  but  sufficient  conditions  and 
proof  have  yet  to  be  found. 

Other  observations  in  all  numerical  examples  we  have  studied  are 
as  follows: 

*  N 

1.  The  optimum  renewal  cycle's  age,  tN  =  2  T.(N)  is  an  increasing 

function  of  N  (Figure  3.6).  *-1 

2.  The  optimum  time  periods  are  such  that 

Tj(N)  £  T.(N+1)  ;  i  =  1,2, ...  ,N  (3.92) 


One  sufficient  condition  that  certainly  must  be  considered  is  the 
fact  that  the  major  repair  cost,  CQ,  is  less  than  the  replacement  cost 


6.0 


x  x  x  x 


I  23456789  10  N 


Figure  3.5-  The  behavior  of  the  minimum  cost 
function  os  N  is  incremented 
(Policy  I-A) 


-90- 


In  order  to  reduce  computations,  initial  values  for  Ti,T2,  • . . 
seem  to  be  helpful. 

This  is  achieved  by  using  Theorem  3.2.  For  instance,  in  the 
Weibull  case,  the  relations  (3.62)  can  be  written  as 


r<TiH)ei  *  r<Ti)ei-l 


Setting  r(t)  =  \Qt  ,  we  get 


i  =  1,2 . N-l  (3.93) 


Vwei 


or 


i  =  1,2 . N-l  (3.94) 


where  0j  is  given  by  (3.75). 

As  an  example,  consider  the  numerical  example  3.1.  The  optimal 
solutions  for  N=1  and  2  are  shown  in  Table  3.1. 


N  C(N,T )  Ti(N)  T2(N) 

1  3.00  10.00 

2  2.96  9.48  0.97 

Table  3.1 

Initially,  we  used  Tt(2)  =  10.0  and  from  relation  (3.94),  a  reasonable 
estimate  is  T2(2)  =  0.91,  which  seems  to  be  close  to  9.48  and  0.97 
respectively . 

Finally,  two  reformulations  of  the  age-dependent  policy  II  are 
presented  without  any  analysis.  Despite  their  additional  difficulties 
caused  by  the  extra  number  of  variables  and  constraints,  some  results 
can  be  obtained  more  easily  by  means  of  interpretation  of  the  Lagrange 


AD- A 1 30  076  SCHEDULING  MAINTENANCE  OPERATIONS  WHICH  CAUSE  OL/ 

AGE-DEPENDENT  FAILURE  RATE . . (U)  POLYTECHNIC  INST  OF  NEW  ' 
YORK  BROOKLYN  DEPT  OF  ELECTRICAL  ENGI . . 

UNCLASSIFIED  B  EBRAHIMIAN  ET  AL.  01  UUN  83  F/G  5/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS-  1963-A 
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multipliers.  Our  objective  is  to  show  that  they  do  exist,  but  no  attempt 
will  be  made  to  derive  equations  satisfying  optimal  solutions. 

The  first  formulation  introduces  the  age  factor  at  the  time  of 
repair  as  a  constraint  by  letting 

°M0i-l  =  Ai-i  •  *  =  2'3 . N  (3-95) 

Thus,  the  expected  cost  per  unit  time  (3.76)  can  be  written  as 

C(N,T)  =  - jP -  (3.96) 


This  cost  function  is  to  be  minimized  in  the  following  nonlinear 
programming  problem 


Find: 

*  *  *  *  * 

Ta . TN'Ai'A2 . AN-1 

To  minimize: 

C(N,T) 

(3.96) 

Subject  to: 

CMei-i  =  Ai-i  i  =  2,3, . . .  ,N 

T  SO 

(3.97) 

N  >1 

N  integer 

The  Lagrange  multiplier  can  be  used  to  analyze  this  problem. 
Mathematically,  this  is  equivalent  to  the  one  we  discussed  previously. 
Lagrange  multipliers  associated  with  constraint  (3.97)  are  found  to  be 


Mi  = 


Ti 

/01r(t)dt 

N 

•  Z  Ti 
j=l  ] 


i  =  1,2, . . .  ,N  (3.98) 
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The  last  formulation  deals  with  the  following  ratio 


Pi  = 


T. 


l 


N 

I  T. 
j=l  ] 


N 


0.  can  be  expressed  as 
8i  *  1  ♦ 


The  long-term  average  cost  per  unit  time  is  given  by 


Cr+(N"1)C0+Cm  Itt*  €TN/iliPm)J0f,itNr(t)dt 
C(N,T)  =  - r - - - 


and  the  resulting  nonlinear  programming  problem  is 
Find .  N  ,Piip2i  ■  •  • 

To  minimize:  C(N,T) 

N 

Subject  to:  I  0.  =  1 

i=l  1 

T  SO 

N  SI 

N  integer 


3.11  MODEL  B 

The  age  factor  for  this  model  may  be  defined  as 


=  €  l  T.  i  =  1,2 . N;  S  0 

1  j=l  J 


(3.99) 

(3.100) 

(3.101) 

(3.101) 

(3.102) 

(3.103) 

(3.104) 
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and  the  failure  rate  for  each  time  period  is 


i-1 

r^tie.,)  =  r(t)  +  €  I  T. 
111  j=l  J 


i  =  1,2, ...  ,N 
t,  Z  0 


(3.105) 


where  r1(t|6i  l)  =  r(t). 

Figure  3.2  is  an  illustration  of  this  case  when  failure  rates  are 
linearly  increasing.  The  long-term  average  cost  per  unit  time  is  found 
by  combining  (3.9)  and  (3.105). 


N  T. 


N  i-1 


c(n,T)  = 


Cr+(N-1)C0+Cm  Z=i/0  V(t)dt+  6Cm  ia.  1  T.) 


N 

■Z  Ti 

j=l  J 


(3.106) 


Taking  its  first  partial  derivatives  and  setting  them  to  zero,  we  get: 


N 


N 


NT.  N  j-1 

[r(T.)  +6ZTJZT,  -  Z  L]r(t)dt  -  61  (T.  1  T) 
1  m?fi  m  j=l  J  j=l  0  j=l  Jm=2m 


Cr+(N-1)Cc 

°M 


i  =  1,2,. ..,N 


(3.107) 


From  which  the  minimum  expected  cost  per  unit  time  is  obtained  by 


*  N 

C(N,T  )  =  Cj^lror.)  +  €1  Tm]  i  =  1,2 . N  (3.108) 

m# 


which  yields 


r(T.)  + 


N 

6  I  T 
m?i 


m 


or 


N 


r(T.)  +  6  Z  T 
J  m^j 


m 


r(T.)  -  t(T.)  =  €  (T.-Tj) 


(3.109) 


(3.110) 
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One  obvious  solution  is  Tj  =  Tm  -  T^.  To  see  if  this  produces  a 
minimum  point,  we  evaluate  |Ey|,  which  has  been  defined  by  (3.61)  as 


icy  =  inwi 

The  second  partial  derivatives  of  R(N,T)  are 


li;  =  W 

=  V 


and  consequently 


J  =  i 


j  t  i 


'd/i  r  si 


r'(Tx)  €  £  ...  6 

6  r'(T2)  6  ...  6 
£  6  r'(T3)  ...  £ 


€  €  €  r'(T.) 


i  =  1,2 . N 


Setting  Tj  =  Tj  =  TN  for  all  i  and  j  yields 


'V  =  si 


r’(TN)  £  6 

6  r’(T*) 


€ 

€ 


6  €  r'(TN)  ...  £ 


€  €  6 
i  =  1,2, ... ,N 


r(TN> 


which  gives 


IPfl  =  r’(TN) 


(3.111) 


(3.112) 


(3.113) 


r 


IDil  =  [r'(TN)}2-62 

=  [r'(T*>  -  €Hr'(T*)  +  £] 


m\  =  [r'(TN)]3  -  3£2r'(TN)  +  2  €3 
=  [r'(T*)  -  €]2lf'(T*)  +  2  6] 


and  in  general 


|Dj|  =  [r'(TN) 


el^Ir'CT*)  +  (i-l)6] 


(3.114) 


In  order  to  have  a  minimum,  it  is  required  that  |D'|  >0  for  all  i 

* 

and  hence  the  sufficient  condition  that  T.  =  Tj  =  to  be  a  minimum  is 
expressed  by 


r’(TN)  >  6 


(3.115) 


Setting  T.  =  TN  in  relation  (3.107),  we  get  an  equation  that  can 


be  solved  for  TN. 


[r(T*)  +  6(N-1)T*]NT*  -  N/0Nr(t)dt 


6(N-1)N  Tl  _  CR+(N'1)Co 
N  '  °M 


or  gN(T*)  s  T^r(T* )  -  r(t)dt  T* 


CR+<(,-1>Co  a 

—rcr. —  -  h(N> 


(3.116) 


The  following  theorem  is  useful  for  finding  the  sufficient  condition 

* 

for  minimum  and  an  estimate  of  N  . 
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Theorem  3.4 

If  (i)  r(t)  is  strictly  increasing  in  N  and  (ii)  >  CQ,  then  the 
solution  to  equation  (3.116)  is  decreasing  in  N  and  unique. 

Proof.  For  the  right-hand  side  of  (3.116),  we  have 

-C  +C 

Ah(N)  =  h(N+l)  -  h(N)  =  (3.117) 

From  hypothesis  (ii),  we  conclude  that 

Ah(N)  *  0  (3.118) 

For  the  right-hand  side  of  (3.116),  we  can  write 
AgN(x)  =  gN+1(x)  -  gN(x)  =  \  x2  £  o 


or  gN+1(x)  *  gN(x)  (3.119) 

which  implies  that  gN  is  increasing  in  N.  From  (3.118),  we  have 
Ah(N)  =  gN+1(T*+1)  -  gN(T*)  <  0 

or  <  gN(TN)  (3‘120) 

By  using  the  implication  of  (3.119),  the  left-hand  side  of  (3.120)  can 
be  decreased  to  get 

*  * 


9N(TN+1'  <  gN(TN' 
and  consequently 


*  * 
tnh  <  tn 

*  * 

Tn  is  unique  since  h(N)  is  independent  of  Tv  and 


■N 


(3.121) 


dgN(TN)  *  *  * 

-NN  =  V(tn)  ♦  6(N-1)Tn  S  0 

dTN 


(3.122) 


Hence,  is  unique  and  decreasing  in  N. 


One  application  of  the  preceding  theorem  is  to  modify  sufficient 
condition  (3.115).  Since  is  the  largest  value  amont  T.'s,  we  have 


r'(Ti)  >  6 


(3.123 


where  Ti  is  the  solution  to  equation  (3.19). 

As  an  example,  consider  a  new  system  having  Weibull  failure  time 
distribution  given  by 


f(t)  =  Xtaexp(-  ) 


ail;  t  >  0;  \  >  0 


Ti  is  given  by  (3.21)  and  hence 


r(a+DcRii/<«+i) 

*  ' 6 


(3.124 


Therefore,  in  the  Weibull  case,  if  the  age  deterioration  factor,  €, 
is  small  enough,  then  the  optimal  solution  is  unique  and  all  major  repair 
time  intervals  are  equal. 

To  obtain  the  optimal  N,  we  set  T-  =  T  in  relation  (3.106)  to  get 


T  6N(N-1)Cm 

CR+(N-l)C0+NCM/^r(t)dt+  - j ^  T2 


We  may  treat  N  as  a  continuous  variable  and  set  the  derivative  of 
(3.125)  with  respect  to  N  to  zero 


3C(N,T)  _ 
3N 


~CR+Co 


eCM 

+  T  =  0 


which  yields 
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N  = 


1 

~T 


I 


2<CR-Co> 


(3.125) 


Let 


N'  = 


rN* 


2(cr-C0) 


6CM 


(3.125) 


Thus  N*  =  [N']  or  [N']+l  depending  on  the  minimum  cost  obtained  from 
(3.108). 

_  .  *  .  _ * 


C(N,T  )  =  CM(r(Tjyj)  +6(N-1)T^] 


(3.126) 


Since  Tx  is  the  largest  value  among  T.'s,  we  can  define  a  lower 


limit  for  N 


* 

N  £ 


i  fo0 
l£j  €Cm 


(3.127) 


For  the  Weibull  distribution  (3.127),  this  becomes 


* 

N  £ 


e 

_ 1 

l/(a+l)  r 

(a+l)CR 

l 

2(Cr-C0) 


(3.128) 


This  estimate  for  N  can  save  computation  time  required  to  evaluate 
the  optimal  solutions. 

Let  us  evaluate  equation  (3.116)  in  case  the  underlying  distribu¬ 
tion  is  Weibull 

-£r  tn+1  +  tn  =  h(N)  (3.129) 

When  ot=l  and  3  the  analytical  solutions  can  be  obtained.  For  a=l, 


we  get 
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n/2 


Similarly  when  a  =  3 


eatll 

3A 


/ 


e2CN-i>2 

9A* 


4h(N) 
3  A. 


1/2 


(3.130) 


(3.131) 


Numerical  Example  3.2 

For  a  Weibuil  distribution 

a  kta+1 

f(t)  =  Ataexp(-  )  a  >  1;  A  >  0;  t  >  0 

with  parameter  A=l,  a=2,  and  maintenance  costs  CR=15,  CQ=5,  CM=1, 

and  the  deterioration  rate  =0.1.  We  applied  gradient  technique  to 

C(N,T)  and  the  optimum  number  of  major  repairs  per  cycle  (including 

*  * 

replacement)  was  found  to  be  N  =7  with  average  cost  C(7,T  )=5.14. 

The  optimum  time  intervals  are  T.(7)=2.02;  i=l,2,..., 7.  Figure  3.7 

illustrates  the  behavior  of  C(NfT  )  as  N  is  incremented.  It  has  a 

decreasing  and  then  increasing  behavior  and  its  minimum  is  unique. 

* 

The  lower  bound  using  (3.128)  yields  N  £6. 

Thus  far  we  have  assumed  that  T.=Tj=TN.  In  our  numerical 
examples,  this  seems  to  be  the  only  optimal  solution.  In  the  Weibuil 
case  with  or=l,  this  is  the  only  solution.  For  a^b,  we  can  see  easily 
that 


C(N,T) |T.=a,  t.=b  =  C(N,T)|T.=b,  T.=a  (3.132) 

which  implies  that  equation  (3.107)  can  have  N!+l  solutions.  If  C(N,T) 
has  a  unique  minimum,  then  the  only  solution  is  T^=Tj;  tfj. 


f 

t 

i 


3.12  AVAILABILITY 


In  our  discussion  on  type  II  policy#  we  assumed  that  replacement, 
minimal  repair,  and  major  repair  times  possess  zero  duration.  However, 
by  interpreting  CR,  C^,  and  CQ  as  constant  times  for  replacement, 
minimal  repair  and  major  repair  costs,  we  can  formulate  the  type  II 
policy  as  steady  state  availability 


A(N,T) 


N 

•  2  Ti 

J=L  J 


N 


N 


\ 


I  T^HN-DC^  If0 


=  (1  *  C(N.T)]'1 


(3.133. 


It  is  seen  easily  that,  from  the  optimization  point  of  view,  maxi¬ 
mizing  the  availability  is  equivalent  to  minimizing  the  long-term  average 


CHAPTER  4 


CONCLUSIONS 


4.1  INTRODUCTION 

We  have  studied  the  optimization  of  two  major  preventive  mainte¬ 
nance  policies  for  repairable  stochastically  failing  systems  with  lifetime 
distribution  function  FjCtlQ.^)  and  failure  rate  r-CtlQ.^),  where  i 
indicates  the  number  of  performed  repairs,  and  the  system's  age  at  the 
time  of  the  (i-l)st  repair  is  represented  by  t.  Our  objective  for 
each  strategy  is  to  optimize  a  set  of  successive  maintenance  intervals 
T1/T2,...,Tn  and  the  number  of  maintenance  intervals  N  for  each 
renewal  cycle  so  that  the  long-term  average  cost  per  unit  time  is  min¬ 
imized.  Each  preventive  maintenance  policy  has  been  classified  into  two 
different  classes,  namely,  state-  and  age-dependent  models.  State- 
dependent  refers  to  the  case  when  the  failure  rate  at  the  time  of  the 
i**1  corrective  maintenance  is  a  function  of  previous  number  of  repairs 
cr  rA( 1 1  ej_ i )  =  r-(t).  In  the  age-dependent  model,  the  system's  age  at 
the  time  of  the  i^1  repair  will  determine  the  (i+l)st  failure  rate.  This 
is  usually  shown  by  r.Ctie^)  =  6._jr(t)  where  0.^  =  1  +  €[age(i-l)]. 
6  is  a  constant  positive  parameter  and  is  called  age  deterioration 
factor.  T^e  underlying  lifetime  distribution  is  Weibull  with  a  strictly 
increasing  failure  rate. 

Our  goai  has  been  efficiency,  and  in  this  context,  many  results 
have  been  obtained,  among  which  only  a  few  are  to  be  mentioned. 
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4.2  POLICY  I 

For  state-dependent  policy  I,  it  was  found  that  the  computational 

algorithm  developed  by  Nguyen  and  Murthy  (1981),  can  be  improved. 

# 

This  is  done  by  finding  a  heuristic  procedure.  We  can  compute  N'  by 

*  * 

which  C(N,®)  is  minimized,  where  N'  provides  a  lower  bound  for  N  . 

This  procedure  is  also  applicable  for  an  age-dependent  in  which  the 
underlying  lifetime  distribution  has  the  Weibull  form. 

The  age-dependent  policy  I  was  optimized  by  simulation  and 
approximation  techniques.  In  both  cases,  using  Weibull  examples,  it  is 
observed  that: 

*  _* 

1.  There  is  a  unique  N  that  minimizes  the  cost  function  C(N,T  ). 

2.  The  optimum  planned  repair  times  are  ordered  as 

Ti(N)  >  T2(N)  >  ...  >  Tn(N) 

and 

T i(l)  >  T i(2)  >  ...  >  TX(N*)  >  TX(N*+1)  >  ... 

Where  in  the  state-dependent  model,  developed  by  Nguyen  and  Murthy 
(1981),  Tx(i);  i  =  1,2,..., N  are  ordered  as 

Ti(l)  >  Tx(2)  >  ...  >  TX(N*)  <  T1(N*+1)  <  ... 

3.  The  results  for  simulation  and  approximation  techniques  seem 
to  be  similar  but  not  identical. 

4.3  POLICY  II 

For  a  general  type  II  policy  and  a  fixed  N,  it  has  been  shown  that 
the  cost  function  C(N,T)  has  a  relative  minimum  and  also  that  the 
optimal  T.  satisfies 
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ri(Ti|0i-l)  <  rj(Tj,ej-l)  ;  1  <  j 

These  properties  are  useful  when  a  optimum  search  technique  is  used. 

The  following  results  are  valid  only  for  Weibull  lifetime  distribu¬ 
tions  and  A-type  age-dependent  policy  II. 

1.  It  is  proved  that  the  optimum  major  repair  times  are  ordered 
as 

TX(N)  >  T(N)  >  ...  >  Tn(N) 

£ 

2.  Numerical  examples  have  shown  that  the  minimal  cost,  C(N,T  ), 

*  * 
as  a  function  of  N,  is  decreasing  for  N  <  N  and  increasing  for  N  >  N  . 

N  * 

3.  IT.,  the  total  duration  of  an  optimal  renewal  cycle,  is  an 
i=l  1 

increasing  function  of  N. 

For  B-type  age-dependent  policy  II,  the  following  results  seem  to 
be  important: 

%  £ 

1.  One  solution  is  T.  =  Tj  =  TN  for  all  i  and  j.  This  is  proved 
sufficient  when  the  age  deterioration  factor,  ,  is  not  large  and  the 
underlying  failure  distribution  is  Weibull. 

2.  If  the  replacement  cost  is  more  than  the  major  repair  cost, 
* 

then  Tn,  the  optimal  repair  interval,  is  decreasing  in  N. 

3.  Our  numerical  example  shows  that  the  minimal  cost  C(N,T  )  is 
convex. 

* 

4.  There  exists  a  lower  bound  on  N  ,  which  can  be  used  as  a 
good  estimator. 

5.  Analytical  solutions  are  obtainable  for  some  Weibull  cases. 

6.  If  C(N,T)  has  a  unique  minimum  with  respect  to  T,  then  the 
only  optimal  solution  is  T.  =  T*;  for  all  i,j. 
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4.4  PROBLEMS  FOR  FURTHER  RESEARCH 
A  few  possible  extensions  are  as  follows: 

1.  For  age-dependent  policy  I,  one  may  prove  or  determine  con¬ 
ditions  under  which  the  optimal  expected  long-term  average  cost  per 
unit  time  C(N,T*),  as  a  function  of  N,  is  decreasing  for  N<N  and 

increasing  for  N>N  .  To  do  so,  equations  (2.97),  (2.103),  and  the 

*  N  * 

behavior  of  the  expected  age  of  the  system,  tM  =  £  y. ,  as  N  increases 

w  i=l  1 

seem  to  be  helpful. 

2.  One  may  consider  extensions  similar  to  (1),  but  using  equations 

(3.51)  and  (3.52)  for  policy  II  (type  A).  In  this  case,  the  increasing 
*  N 

property  of  t.r  =  £  T.(N),  as  a  function  of  N,  must  be  proved,  for 
N  i=l  1 

type  B  of  policy  II,  equation  (3.116)  is  essential. 

3.  We  have  assumed  that  the  replacement,  major  repair,  minimal 
repair  and  breakdown  costs  are  constant.  One  may  take  a  more  prac¬ 
tical  approach  by  defining  these  maintenance  costs  to  be  functions  of 
elapsed  operating  time  of  the  system.  In  this  case,  it  might  be  easier 
to  minimize  the  expected  total  cost  rather  than  the  expected  long-term 
average  cost  per  unit  time. 

4.  Our  models  can  be  generalized  by  assuming  non-zero  and 
random  down-times  for  replacement  and  repairs. 

5.  Another  generalization  is  to  have  failure  rates  increasing  in 
the  number  of  major  repairs  and  real  age  of  the  system  simultaneously, 
where  we  have  studied  them  separately. 

6.  Policy  II  can  be  modified  by  including  a  constant  or  minimum 
up-time  constraint  where,  in  policy  I,  average  up-time  may  be  restricted. 


7.  It  was  assumed  that  the  minimal  repairs  always  restore  the 
system  to  operation.  One  may  reformulate  the  policy  II  such  that  failures 
are  corrected  by  minimal  repairs  with  probability  p  £  0. 

8.  It  may  be  possible  to  find  parameter  sets  for  which  the  optimal 
policy  permitting  intermediate  repairs  has  no  such  repairs. 


APPENDIX 


MULTIDIMENSIONAL  MINIMIZATION 
BY  CONJUGATE  GRADIENT  METHOD 


Conjugate  gradients  are  among  the  advanced  search  methods  to 
calculate  an  unconstrained  minimum  of  a  real  function  with  several 
variables.  The  underlying  method  is  found  in  an  article  by  R.  Fletcher 
and  C.  M.  Reeves  (1964). 

The  method  is  used  to  find  X  (a  column  vector)  which  minimizes 
C(X)  by  a  sequence  of  moves  from  an  initial  point  XQ  to  a  new  point 
X1(  then  to  X2,  and  so  on,  where  the  gradient  VC(X)  =  G(X)  is 
available  analytically.  Proceeding  from  an  arbitrary  initial  search  point 
XQ,  we  locate  a  sequence  of  points  that  are  successively  closer  to  the 
minimum  by  the  relation 


where  X.  is  the  current  point  and  a.  is  a  positive  scalar  that  defines 
the  distance  between  X.  and  Xi+1  along  the  search  vector  (a  column 
vector).  Notice  that  the  minimum  along  Sj  will  occur  where  the 
gradient  of  C(X.+p  will  be  normal  to  Sj.  Although  the  algorithm  is 
supposed  to  find  the  optimum  of  a  quadratic  form,  it  is  in  fact  appli¬ 
cable  to  a  broader  class  of  functions.  The  reason  is  that  any  function 
C(X)  in  the  neighborhood  of  the  required  minimum  X*,  can  be 
approximated  closely  by  the  first  three  terms  of  its  Taylor  expansion 
and  consequently  by  a  quadratic  form. 
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Given  XQ  and  G(Xq)  ,  the  conjugate  gradient  method  can  be 


described  as 

follows : 

STEP  0: 

Set  i  =  0;  SQ  =  -G(Xq). 

STEP  1: 

Form  C(a.)  =  C(Xj  +  a.S.). 

STEP  2: 

Find  a*;  the  value  of  a.  which  minimizes  C(a.) 

STEP  3: 

Obtain  new  point  X.+1  =  Xj  +  ot.S.. 

STEP  4: 

Determine  new  gradient  G(X.+1). 

STEP  5: 

Evalutae  C(X.+1). 

STEP  6: 

To  test  for  optimality,  STOP  if  either 

a.  ACj  =  C(X.+1)  -  C(X.)  £  0. 

b.  (G(X.+1)TG(Xi+1)  £  6. 

c.  i  ^  n. 

STEP  7: 

<g(x.+1))tg(x.+1) 

Compute  £.  =  — - - 

1  (G(X.))IG(X.)  " 

STEP  8: 

Find  the  new  search  vector  Sj+^  =  -G(X.+1)  + 

STEP  9: 

Set  i  =  i+1. 

STEP  10:  Go  to  STEP  1. 

where  6  >  0  is  a  predetermined  small  number  and  n  is  the  maximum 
allowable  number  of  iterations.  Notice  that  Step  2  usually  requires  a 
one-dimensional  search. 

We  have  used  this  algorithm  to  find  the  optimum  scheduled  repair 
times  for  age-dependent  policy  I  (approximation  and  age-dependent 
policy  II,  where  underlying  lifetime  distributions  had  Weibull  forms. 
Despite  the  existence  of  nonnegativity  contraint  T  £  0,  the  conjugate 
gradient  converges  to  the  minimum  point.  Step  2  was  carried  out  by 
employing  the  Golden  Section  search  technique,  where  the  initial  inter¬ 
val  of  uncertainty  for  a.  was  restricted  to  be  positive. 
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Slightly  less  than  1.5  minutes  of  execution  time  was  needed  to  find 
the  optimum  planned  repair  times  for  approximated  age-dependent 
policy  II  with  N=2.  Using  simulation,  the  same  problem  requires 
slightly  greater  computation  time,  but  when  N  is  large,  simulation 
becomes  more  efficient  than  approximation.  Much  less  execution  time 
was  needed  for  age-dependent  policy  II  (Type  A).  For  a  typical 
example,  less  than  one  minute  of  computer  time  was  spent  to  minimize 
C(N,T)  by  finding  T*  when  N  took  on  values  from  1  to  20. 

Here,  we  list  FORTRAN  subprograms  written  for  approximated 
age-dependent  policy  I  and  type  A  of  age-dependent  policy  II.  It  is 
assumed  that  the  underlying  lifetime  distribution  has  Weibull  form: 

A  6. 

fi(ti0i-i)  =  Vi-ieyp(_  — );  \> >  °' a  * 1  >  °- 

where  0J  is  given  by  relations  (2.69)  and  (3.75)  for  policies  I  and  II 
respectively . 


nn  n  nnnno  n  nnnnnn 
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SUBfiOOTINE  GBADII(M,T,EPSrCBNICO,CMLAH,  ALPHA,CfG,TSUN) 
DIHENSXON  T(N),G(N)  rC0EFT  (40) 


GBADIX  BOOTIES  EVALUATES  THE  GRADIENT  OF  THE  EXPECTED  LONG-TEBH 
A71BAGE  COST  FUNCTION  C(NfT),  BH  EH  THE  AGE-DEP EH  DENT  POLICY  II 
IS  ONDEB  CO US I DEB AXIOM  (SECTION  3)0).  THE  ONDEBLYING  FAILOBE 
El STB IBOTION  HAS  THE  HEIBO FOBH  WITH  ALPHA  AND  LAB  AS  SCALE 
AND  SHAPE  PABAMETEBS  BESPECT IYELY.  THE  AGE  FACTOB  IS  GIVEN  BY 
EQUATION  (3.75). 

V ABIABLES  THAT  ABE  GIVEN  ABE: 

H  :  NOHBEB  OF  HAJOB  BEPAIBS  INCLUDING  BEPLACEMENT  IN  A  CYCLE. 

T  :  A  1-DIM  ABBAY  OF  SIZE  N  BEPBESENTING  THE  SCHEDULED  BEPAIB  TIMES. 
EPS:  AGE  DETEBIOB ATION  FACTOB 

CBNICO  =  CB*  (N-1)  CO  :  BEPLACEMENT  AND  MAJOR  BEPAIB  COSTS  IN  A  CYCLE. 
CHLAH  »  CM  *LAH  :  (MINIMAL  BEPAIB)  ♦  (HEIBULL'S  SCALE  PABAHETEB). 

V ABIABLES  HILL  BE  BETUBNED  ABE  : 

C:  COST  FUNCTION  C(N.T) 

G  :  A  1-DIM  ABBAY  BEPBESENTING  THE  GBADIENT  OF  C(NrT) 


ALPHA  1  »  ALPHA  *1.0 
TSOH  *  0.0 
DO  10  I  ■  1 , N 
TSOM  *  TSUM  ♦  T (I ) 

10  COK1INUE 
D2  *  0-0 
DO  15  I  =  2.  N 
D2  *  B2  ♦  T  (I )  ♦*  ALPH  A 1 
15  CONTINUE 

COEFT(I)  *  1.0 
M  =  H  -  1 
DO  201®  1  ,  M 

COEFT  (I  *■  1)  =  COEFT(I)  *  EPS*T  (I) 

20  CONTINUE 
D1  *  0.0 
DO  25  I  *  1,  N 

01  >  DM  COE  FT  (I)  *T  (I)  ** ALPHA  1 
25  CONTINUE 

B  »  CINICO  ^  CHLAH  *D1/ALPHA1 
B 1  *  CHLAH*  (T  (1)  **ALPHA  ♦  EPS*D2/ALPHA  1) 

C  *  B/TSUH 
G  ( 1)  »  (B  1-  O/TSUH 
DO  30  I  *  1  ,  H 

B2  «  COEFT(IM)*T  (M)  **  ALPHA-COEFT  (I)  *T(I)  **  ALPH  A- IPS*T  ( I*  1)  ** 
* ALPS A 1/ALPHA1 
B2  ■  CHLAH  *  B2 
G(I  ♦  1  )  *  G  (I)  ♦  B  2/TSUH 
30  CONTINUE 
BETUB1 
END 


nnnnnon  n  nnnonnn  nonnon 
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SUBBOGTINE  GRADI  (C) 

COMMCI  /BLT/T  (20)  /BLN/N/BLTH/THETA (  20)  /BLY/Y  (20)/ 

$  BLG/G  (2 0) /B LL/L  AMO/BL A/ ALPHA/BLC/CR , CB/BLC0/C0/BL  A1 /ALPHA  1 / 
S  BIS/EPS 

DIMENSION  7(2  0)  ,FC{ 20)  ,0(20,20)  ,FCT  ( 20)  , B  (20) 

REAL  I  AMO  ,  L  ( 2  0) 

ALB*L AMO/ALPH  A1 


BOOTIME  GBAO  I  E7AL0ATES  THE  GHADIENTOF  THE  EXPECTED  LONG-TEBM  AVERAGE 
COST  C (N, T)  ,  WHEN  APPBCXIBAT  ED  AGE-DEPENDENT  TYPE  II  POLICY  (DISCUSSED 
IM  SECTION  2. 9  )  IS  APPLIED. 

THE  UNDERLYING  FAILOBE  DISTRIBUTION  HAS  THE  SEIBDLL  FORM  GIVEN  BY  EQ. 
2.116  AND  THE  AGE  FACTOR  IN  THIS  CASE  IS  DEFINED  BY  EQ.  2.117 
OTEEB  VARIABLES  ABE  : 

B  :  NUMBER  OF  MAJOR  REPAIRS  INCLUDING  REPLACEMENT  IN  A  CYCLE.  f 

T  :  A  1-DIB  ARRAY  OF  SIZE  N  REPRESENTING  THE  SCHEDULED  REPAIR  TIHES. 

EPS:  AGE  DETEBIOB ATION  FACTOB 
C:  COST  FUNCTION  C(N,T) 

G  :  A  1-DIM  ARRAY  BEPBESENTING  THE  GRADIENT  OF  C(N(T)  , 

ALPHA  1  *  ALPHA*! 

AGE  =  SOM  OF  Y'S 

VABIABLE 

I 
Y 

THETA 
D 

a 

L 
G 


i* ALB**  ( 1.  0/ALPH  A  1) 

XH*N*1  (1) 

Y(1)=IHEA(1,XB)/i 
THETA  (1)  *1.0 

V  ( 1 )  *AREA  (2,XB)/H 
ALE* ALB*EPS 
00  5  1*2, I 
B*l** ALPHA  1  +  ALB*Y  (1-1) 

«*■♦♦  (1.0/ALPHA1) 

XB*1*I  (I) 

Y  (I)*ABEA(1,XB)/B 

THETA  (I)  *TEET A  (1-1)  *EPS*Y  (1-1) 
V(I)*AHEA(2,X8|/{H*THETA{I)  ) 

5  CCNTIBOE 
DO  10  1*1,1 
PC  (X)  *FBAB  (T  ( I)  ,1) 

D(I,I)»FC|I) 

10  COBTIIOE 


EQ.  NO. 

(2. 124) 
(2. 123) 
(2.  1 17) 
(2. 112) 
(2.  1 1S) 
(2. 1 14) 
(2.  1 15) 


J 


-112- 


DO  15  J=1,H1 

ai*o.o 

J1*JO 

DO  15  I*J1,H 
i1*&1*D(I-i,J) 

D(I.J)=-BPS*i1*V(I) 

15  com  i  sob 

DO  20  1*2, S 

FCT (I) *T (I) ** ELPHBl *IC (I) 

20  cosmos 

B(H>  *CB*L4H0*THET&(H)*T(H)  **ILPHA 

a  i  s)=fl(N)  *pc(») 

DO  30  1*1, HI 
£*0.0 
ST*  0.0 
DO  25  H*I,N1 
E*E*D  (9,1) 

ST*ST*E*FCT(M  M) 

25  COHTIBOE 
Sl*SI^iLB 

PT*L4H0*THETA  (1) *1 (I) **HLPH4 
FT*FT*PC (I) 

B(I)=CB*(FT*ST) 

30  CCS1IN0E 

EG E*  (1HET&  (H)  -  1.  0  )/EPS  ♦  I  (R| 

DO  35  1*1, H 
L (I) =0. 0 
DO  35  J=I,H 
L  (I)  =1  (I)  *D(J,I) 

35  COITIKOE 

B1*CB*(H“1)*C0*,N+CB 
R2*0. 0 
DO  40  r*t,M 
B2*B2*D  (1,1) 

40  COHTIBOE 

B1*B1-CB*B2 
C*B 1/iGE 
DO  45  1*1, H 
G(I)*B(I)-L(I)*C 
G(I)*G(I)/iGB 
45  COHTIBOE 
BETOBB 


n  n  n 
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FU  HCTIQN  ABE A  (K ,  XB) 

DIHESSION  B ( 40, 40) 

ABZA(K,X)  PERFORMS  HOMES IC  XL  INTEGBATIOH  IN  ITEBVAL  (0 
EXP(-X**(ALPHAM)  )  IP  K=1 

X**(ALPHAM)*EXP  (-X**  (ALPHAH)  )  IP  K=  2 


DO  131  1=1,40 
DO  13  1  J=  1,40 
B (I, J) =0. 0 
131  CONIIBUE 
XL*0. C 
DEL=XB-XL 

B(1,1)=DEL*(FD(K,X1)  »-FD(K,XB)  )  /2.  0 
B(1,2)=B(1,1)  *DEL*FD  (K,  (XL>XB)/2.  0) 
B(1,2)=B  (1,2)  /2.0 
8(2*1)=4*B  (1 , 2)-B  (1,  1) 

B  (2,  1)  =B  (2, 1)  /3.0 
J=3 

5  D=DEL/  (2**  (J  -  1) ) 

X=XL-D 
NH»2**  (J-2) 

saH=o.o 

DO  10  1=1, HH 

X=X*2*D 

SOH*SOS*FD  (K ,  X) 

10  CO  NT  I  BOB 

B ( 1 , J) *B { 1 , J- 1)/2.0*D*SaM 
Ll  =  1 

DO  15  L=2,  J 

KK»J-IH 
Ll*l 1*4 

B(£,KI)=L1*B(L-1,KKf  1)  -B(L-1,KK) 

B (I, KB)  »B  (L,KK)/ (Li- 1) 

15  COMTXIOE 

B1  =  (B  (J,  1)-B(J-1,  1)  )/B  (J,  1) 

B 1  =  ABS  (B 1 ) 

IF  (B1.LT.1. OE-6)  60  TO  20 

J=J*1 
60  TC  5 
20  ABE1=E  (J,  1) 

BETOBI 

BUD 


X)  FOB 


n  n  o  on 


FOICTIOI  FBA  B  (X  ,  K) 

COH DC I/BLL/L AMQ/BLA 1/ALPHA 1/BLTH/THET A (20 ) 
BEAL  I 4 HO 


FB IB  BILL  BETOBN  THE  VALOE  OF 
IXP(-X**(ALPHA*1)  *LAHO*THETA/(ALPHAM)  ) 

X1*LAtO*THETA  (K)  /ALPHA  1 

X2*X** ALPHA  1 

X3*X1*X2 

IF  (X3  .LT.  1-0-75)  GO  TO  1 
FB ABSEXP (-X3) 

BET05I 
1  PB  A  E  ®  1 .  0 
BETOBI 
BHD 


FOICTIOI  FD (Kf X) 
COHRCB  /BL  At  /  ALPH  A 1 


FD(K.X)  BILL  BETDHM  THE  FALOE  OF 
EXP (-X** { ALPHA ♦ 1 )  ) 

X**  (ALPHA*- 1)  *EXP  (-X**  (ALPHA  ►  1)  ) 


IF  (X  -LE.  1.0E-75  .  AIO.  K  -EQ.  2)  GO  TO  4 
IF  (X  -LB-  1.  OB-75  -  AID.  K  -EQ.  1)  GO  TO  5 
X1»X*«ALPHA1 

IF  (XI  -GT-  1.0E-75)  GO  TO  3 
X2* 1- 0 
GO  TO  2 

3  IF  (XI  -GT-  170-0)  GO  TO  4 
X2-1XP(-X1) 

2  IF  (K.EQ.  1 )  GO  TO  1 
FD-XMI2 
BETOBI 
1  FD-X2 
BETOBI 

4  FD»  0.0 
BETOBI 

5  FD-1-0 
BETOBI 
BID 


IF  K=  1 
IF  K=  2 
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